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ABSTRACT 
Title of Thesis: Predicting the ignition and burning rate of wood in the 
Cone Calorimeter using an integral model. 


Name of degree candidate: Michael John Spearpoint 


Degree and Year: Master of Science in Fire Protection Engineering, 1999. 
Thesis directed by: Dr. James G. Quintiere, Professor, Department of Fire - 
Protection Engineering 


This study compares ignition and burning rate measurements of wood in the Cone 
Calorimeter with a one-dimensional integral model that describes the transient pyrolysis 


of a semi-infinite charring solid subject to a constant radiant heat flux. 


Four species of wood were exposed to a range of incident heat fluxes both with their 
grain parallel and perpendicular to the incident heat flux. The time to ignition and 
burning rate measurements obtained from the Cone Calorimeter were used to derive 
characteristic properties of the materials. These properties were used as input to the 


integral model to compare its predictions with the experimental data. 


Thermocouples were embedded as various depths in the samples. Temperature 
measurements were used to determine the progress of the thermal penetration wave and 
char depth through the material. The measurements were compared to predictions from 


the integral model. 
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NOMENCLATURE 


a pre-exponential factor for char fraction relationship [-] 
a thermal diffusivity, [m’/s], absorptivity [-] 
A area, [m’] 


B Ratio of convective gain and radiative loss with incident heat flux, [-] 


io) 


specific heat [J/kg.K] 
ignition constant, [-] 


G 
D diameter (of cylindrical flame), [m] 


ro) depth [m] 

A dimensionless depth, [-] 

é flame emissivity, [-] 

¢, @' char fraction, [-] 

f grain orientation coefficient, [-] 


AH. heat of combustion [J/kg] 

AH, heat of vaporisation or heat of pyrolysis [J/kg] 
n height, [m] 

h heat transfer coefficient [W/m’.K], enthalpy [J] 


thermal inertia, kpc, [J?.m~*.K*.s"] 


~~ 


k thermal conductivity, [W/m.K] 
K absorption coefficient, [m7] 
i heat of gasification, [J/kg] 
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Q 


mean beam length of flame, [m] 

mass, [kg] 

dimensionless mass loss rate, [-] 

power factor for char fraction relationship [-] 
pressure, [Pa] 

pi, [3.141] 

peak mass loss rate ratio [-] 

power output, [W] 

heat flux, [W/m7] 

density, [kg/m*] 

specific gravity, [-] 

temperature, [°C] or [K] 

time, [s] 

dimensionless temperature, [-] 
dimensionless time, [-] 
Stefan-Boltzmann constant [W/m?.K] 
shrinkage factor, [%] 

moisture content, [%], internal energy [J] 
velocity, [m/s] 


volume, [m°] 


constant used in pure convective loss ignition analysis, [-] 
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final 
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gas 

incident 
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measured 
constant pressure 
peak 

radiation 
surface, steady (when subscripted to 6) 
steady-state 
thermocouple 
vaporisation, constant volume 
virgin wood 
across grain 


along grain 


XXXVIi 


tr transitional 


Superscripts 
( )” _ perunitarea 


(-) per unit time 


XXXVIli 


aiserith 


initial, eesti 


boi dhe 


iTri¢e 
” » . 


p Aad 


‘qwoontabep, omelet wolaviae 
Vir git wr Ther | 


OCruULLi ym 


rote grean 


I rt 
iv F e¢ 
i Poy oN “ee ; 
. ¥ 
¥ ye 7 
: Dt ie be 7 as, 
7 . ~~ ae 
t s & pice ar 


1. INTRODUCTION 

“We are not beginners 

we will not be fooled 

the times have been our teachers & they 

teach an iron rule, we do our 

best work in the dark, we do our 

thinking on the run” 

A Fire is Burning, Oysterband (Holy Bandits, 1993) 

1.1 Background 
Ever since prehistoric times humans have known that wood burns. The ability of wood to 
burn has been both a benefit and a problem for humans. The capability to predict the 
burning rate of wood in modern times has become increasingly important as fire safety 
engineering moves toward a performance-based approach to building design. Computer 


based fire and hazard models require the burning rate of materials to be specified as 


input. 


The pyrolysis behaviour of solid materials can be divided into two types: non-charring 
and charring. Non-charring materials burn away completely leaving no residue and can 
be modelled using theory similar to flammable liquids. In contrast, charring materials 
leave relatively significant amounts of residue when they burn. The pyrolysis of charring 
materials such as wood is a complex interplay of chemistry, heat and mass transfer. 
Charring materials must be modelled in terms of a pyrolysis front penetrating into the 


material with an increasing surface temperature and without a well-defined steady state. 


The purpose of this study is to examine a one-dimensional integral model for charring 
materials (in this case wood) by comparing the model with data obtained from an ignition 
and burning rate apparatus (the Cone Calorimeter). This study will show how we might 
predict the ignition and burning characteristics of a charring material using a limited 
number of experimental measurements to obtain properties that are required by the 
model. Although some material properties can be obtained from the literature, several 
properties are derived from the time to ignition measurements and through an iterative 
process of comparing the model predictions from the integral model with experimental 


test data. This study demonstrates a methodology for obtaining these properties. 


1.2 Structure and thermal decomposition of wood 

1.2.1 Growth and structure 

There are many species of wood and they are used by society for a wide range of - 
construction, utensil and decorative purposes. Woods are subdivided into two botanical 
categories: softwoods and hardwoods and the physical structure as well as the chemical 


composition is considerably different between these two. 


Janssens [1] in his thesis gives an excellent overview of the composition of wood and its 
thermal decomposition. Much of the description given is this section is taken from his 


work supplemented by information from Drysdale [2] and Cholin [3]. 


The macroscopic structure of softwood and hardwood stems is very similar. The stems 


consist of a core of wood (xylem) covered by a protective layer of bark. The xylem 


consists of long fibre-like cells that are oriented preferentially in one direction (vertical) 
referred to as the grain. Water, minerals and nutrients are taken up by the roots and are 
transported through the outer part of the xylem to the leaves. Photosynthesis in the leaves 
uses the water and CQ> (from the atmosphere) to form various sugars. A solution of the 
sugars into water (also referred to as sap) is moved to various parts of the tree through the 
inner layer of the bark (phloem). The growth process is limited to a thin layer between the 
xylem and the phloem, called the cambium. New cells are created on either side of the 
cambium through division of existing cells. Thus, both the phloem and the xylem are 


growing. The cambium forms a ring with a continuously increasing diameter. 


After a number of years, the cells in the inner part of the xylem die. The wood in this part 
of the xylem is referred to as heartwood. The remaining (outer) part of the xylem is called 
sapwood. From the moment the first heartwood is formed, it expands together with the 
cambial layer. Therefore, commercial lumber usually consists of a mixture of primarily 
heartwood and some sapwood. The mechanical properties of heartwood and sapwood of 
the same stem are very similar. However, heartwood is often denser than sapwood. In 
addition, heartwood contains a variety of products of sugar decomposition originating in 
the cambium. These products can be removed by physical or chemical extraction 
processes and are therefore commonly referred to as extractives. The extractives clog tiny 
passages in the cell walls. Consequently they significantly reduce the permeability of 
wood. The chemical composition of some extractives is very different from that of the 
main components of wood; cellulose, hemicellulose and lignin. Thus, because of the 


aforementioned physical and chemical phenomena associated with the presence of 


extractives, the fire behaviour of heartwood may be somewhat different from that of 


sapwood. 


The elemental composition of dry wood consists of nett 50% carbon, 6% hydrogen, 
44% oxygen and small amounts of nitrogen and some inorganic compounds. The 
principal elemental constituents are combined into a number of natural polymers; 
cellulose, hemicellulose and lignin. Cellulose consists of a large number of glucose 
molecules joined together in a chain-like polymer. The molecules align themselves into 
bundles (microfibrils) which provide structural strength and rigidity to the cell wall. 
Glucose, (C¢H)20¢), 1s the principal sugar generated by photosynthesis. The chemical 
formula of cellulose is (C6H}09Os), i.e. one molecule of water is generated for every pair 
of glucose molecules linked together. The degree of polymerisation, m, may be as high as 
30,000. Various other sugars produced in the leaves are combined to branched-chain 
polymers called hemicellulose. The degree of polymerisation is generally only a few 
hundred. Lignin is a stable high molecular weight polymer that is phenolic in nature. It 
acts as a binding agent within and between cell walls. The lignin content, in general, is 


significantly higher in softwoods as illustrated in Table 1. 


Type Cellulose Hemicellulose Lignin 
Hardwood | 40-44 23-40 18-25 
Softwood | 40-44 20-32 25-35 


Table 1. Chemical composition of dry wood in percent mass. | 


Wood normally contains a quantity of moisture present by different physical 


mechanisms. Much of the space in a live tree is occupied by water, called “free water”. 


When a tree is cut, the free water is gradually replaced by air. Wood also contains water 
that is hygroscopically bound to the cellulose molecules through a weak chemical bond. 
This “bound water” will leave the wood at a rate determined by the temperature and 
relative humidity. The moisture content of wood below the range of 5 to 7 per cent is 
called the “moisture of constitution”. This water forms the bonds between the adjacent 
cellulose molecules and cannot be removed without the physical destruction of the wood. 


Thus, even “dry” wood always contains some moisture. 


Due to its grain, wood is an anisotropic material i.e. it does not have the same structure in 
its orthogonal directions. Thus the structure of wood can be considered parallel to the 
grain, tangential to the grain and radial to the grain. As will be discussed later in more 


detail, the properties and burning characteristics vary with the grain orientation. 


1.2.2 Burning characteristics 
As a result of the complex chemical and physical structure of wood, its ignition and 
burning characteristics are also complex. Consequently, only a simplified description of 


these mechanisms is given here. 


Consider a slab of wood with one surface uniformly exposed to a constant external 
incident heat flux and a source of ignition (a pilot) close by. The homogeneity of the 
incident heat flux to the exposed surface allows the heat and mass transfer through the 
slab to be considered one-dimensional. The incident flux pyrolyses the surface of the 


wood thus releasing volatile fuel gases. The mass flow rate of the volatiles is a function 


of several factors including the intensity of the incident energy and the orientation of the 
grain. If the mass flow rate of the volatiles is sufficient to reach the lower flammable 
limit of the fuel/air mixture then ignition occurs. At this instant the heat flux to the 
surface of the wood is a combination of the external flux plus the flux from the flame. 
The rate of heat release rapidly rises to a maximum then a char layer gradually builds up 
as the pyrolysis front moves inward. The char layer forms an increasing thermal 
resistance between the exposed surface and the pyrolysis front resulting in a continuously 
decreasing rate of heat release after the first peak. Experimental data are available 
showing that if the slab is sufficiently thick, its rate of heat release eventually reaches a 


more or less steady value. 


The char layer begins to break down rapidly at reiterate above 300 °C. The char 
layer also shrinks and pressure gradients (i.e. stresses) are set-up within the material. 
Small cracks appear on the surface, perpendicular to the grain, and these cracks allow 
volatiles to escape more easily. The cracks gradually widen as the char layer deepens 


leading to the characteristic ‘alligatoring’ patterns that are frequently observed. 


After the volatiles have been exhausted, flaming ceases and a solid char residue remains. 
The char continues to burn in a smouldering mode. Prior to that, char oxidation is usually 
minimal since the flame prevents diffusion of O to the surface. Therefore, heat release 
rate and related quantities measured during the flaming phase are predominantly those of 
the volatiles. The main constituent of wood char is carbon, so that its net heat of 


combustion is around 32 MJ/kg. It has been found for various species of dry wood that 


the char yield is 33% + 3% of the original mass of wood. With an average net heat of 
combustion of wood of 17 MJ/kg, the mean heat of combustion measured in the Cone 
Calorimeter during the flaming mode should then be about 10 MJ/kg. The values 

(§ 4.4.1) actually measured are somewhat different due to variations in the char yield and 


moisture in the conditioned samples. 


The three main components of wood have quite different thermal degradation 
characteristics. This is illustrated by thermogravimetric analysis, showing that the 


constituents decompose to release volatiles over different temperature ranges [2], 


typically: 
Cellulose 240-350 °C 
Hemicellulose 200-260 °C 
Lignin 280-500 °C 


Consequently, the thermal degradation characteristics of wood shift towards higher 
temperatures with increasing lignin content. This explains why the surface temperature at 
ignition is significantly higher for softwoods than for hardwoods (see § 3.2.3). In 
addition, at temperatures in excess of 400-450 °C lignin decomposes to volatiles for 
about 50% of its mass and is therefore responsible for most of the char. In contrast, 
a-cellulose leaves only 5% char after prolonged heating at 300 °C. The charring rate of 
wood is also very sensitive to the presence of inorganic impurities, such as fire retardants, 


because they affect the chemical kinetics of the pyrolysis process. 


The heat of combustion of wood measured in an oxygen bomb, which includes the heat 
released by the char, increases with lignin content. This may not be relevant to real fire | 
situations of interest for this work, since the heat released under such conditions is 
primarily that of the volatiles. However, the net heat of combustion of the volatiles is also 
slightly higher for lignin than for the other constituents. Parker [4] measured 14.7 MJ/kg 
for lignin, 13.8 MJ/kg for cellulose and even lower values for some hemicelluloses of 
Douglas fir. Still, since more energy is required to generate the volatiles, it is expected 


that fire performance of wood improves with increasing lignin content. 


As a result of the grain structure of wood, not only do the material properties vary but 
also the mechanisms by which the volatiles are able to exit the surface of the wood and 
thus the ignition and burning characteristics. Consider the xylem cells as long “tubes” and 


two scenarios where the incident heat flux is either parallel or perpendicular to these 


tubes (Figure 1) 


Incident Incident 
heat heat 
Preferred direction of 


travel of volatiles 


(a) Incident flux parallel to grain. (b) Incident flux perpendicular to grain. 


Figure 1. Incident flux and grain orientation scenarios. 


Volatiles generated just below the surface of the unaffected wood can travel more easily 
along the grain compared to at right-angles. The appearance of jets of volatiles and 
flames at the ends of a burning log is evidence for this. In case (b), the cell wall impede 
the flow of volatiles to the exposed surface and thus the walls have to decompose to 
allow a sufficient mass flux of fuel to achieve the lower flammable limit. This 
decomposition process requires additional energy and thus we might expect (and this 
study demonstrates, at least for low heat fluxes) that the ignition of wood is more easily 


achieved when the wood is exposed with its grain end-on 1.e. case (a). 


1.3. The Cone Calorimeter 

1.3.1 | Background 

The Cone Calorimeter [5] is a standard apparatus for measuring the ignition and burning 
characteristics of materials. The apparatus was developed during the 1980's at the 
National Institute of Standards and Technology (NIST) formerly the National Bureau of 
Standards (NBS). The Cone Calorimeter at the University of Maryland was installed and 
tested by Woodford [6]. It is not the intention of this study to provide ceils of the theory 
and operation of the apparatus as this can be readily found in the literature. For example, 
the construction and operation of the apparatus is detailed in NFPA 264 [7] and ASTM 
1354 [8]. The University of Maryland Cone Calorimeter and its accompanying software 
was designed and built during or prior to 1988 and thus meets earlier versions of the 
relevant standards than those currently in use. Figure 2 shows a schematic diagram of the 
Cone Calorimeter taken from reference 8. The University of Maryland Cone Calorimeter 


includes the chamber. 
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NOTE 1—All dimensions are in milimetres. 
Note 2—° indicates a critical dimension. 


Figure 2. Schematic of the Cone Calorimeter (from ASTM 1354, 1995). 


Samples are mounted in a retainer frame. A truncated cone shaped electric heater is used 


to impose an external incident heat flux onto the surface of the sample. An electric spark 
igniter can be used to initiate piloted ignition. The time to ignition ¢,, is recorded by the 


operator. The apparatus uses a load cell to measure mass loss rate per unit area m” . 


Oxygen consumption calorimetry is used to obtain the rate of energy release per unit area 


Q” . The combustion products from the burning sample are collected in the hood and 
passed along the duct. The relative decrease in the oxygen concentration in the 
combustion gases compared with normal ambient is recorded. The mass flow of the 
combustion gases is obtained from the measurements of the temperature and pressure 
difference across an orifice plate located in the stack. A thermocouple is used to measure 
the gas temperature and a pressure transducer is used to obtain the pressure difference. 


Thus, the amount of oxygen consumed per unit mass of combustion gas can be found. By 
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using the fact that the energy release of a burning material per unit mass of oxygen 


consumed is almost a constant for most materials [9], the energy release rate can be 


obtained. The instantaneous heat of combustion AH, can be obtained from 


Equation 1 


and it has been observed that the heat of combustion is generally constant for a material 
undergoing flaming combustion. In addition, the smoke density across the duct is 
determined by using a receiver to measure the obscuration of a laser as smoke passes 


through its optical path. 


The mass loss rate, rate of heat release (rate of energy release or also referred to as the 
burning rate), instantaneous heat of combustion, smoke obscuration and other data were 
all automatically recorded as functions of time by a commercial data acquisition system 
connected to the Cone Cree For the majority of the tests, the scanning rate of the 
data acquisition system was 5 s. However, smaller scan rates of 2 s and 1 s were used in a 
limited number of the experimental tests. Note that the current (1995) version of the 
NFPA standard [7] requires a scan rate of 2 seconds or less as opposed to 5 s or less in 
earlier versions. In addition, the definition of sustained flaming is 4 seconds rather than 


the earlier definition of 10 seconds which is used in this study. 


The data collected by the Cone Calorimeter data acquisition system was post-processed 
using a program provided with the system. Finally, this processed data was imported into 
a commercial spreadsheet package for analysis. 
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1.3.2 Calibration and set-up 

In preparation for the study, a retainer frame was modified to allow thermocouples to be 
located inside the test sample. The heat flux from the cone heater was calibrated against a 
calibrated total heat flux gauge. A second standard reference total heat flux gauge was 
used to verify the calibration. Thus the temperature of the cone heater was selected to 


provide a known incident heat flux (Figure 3). 
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Incident heat flux, q”, [kW/m?’] 


700 800 900 1,000 


Cone heater temperature [°C] 


Figure 3. Comparison between the temperature setting on the Cone Calorimeter and the 
calibrated total heat flux gauge. 


Prior to each day of testing the Cone Calorimeter was calibrated according to the 
procedure laid down by Woodford [6]. The load cell was calibrated between 0 g and 


250 g in increments of 50 g using standard weights. The pressure transducer in the duct 
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was zeroed. The oxygen analyser was calibrated with nitrogen (i.e. 0% oxygen) and 
ambient (i.e. 20.95 % oxygen). The smoke extinction laser system was calibrated using 
0.3 OD/m and 0.8 OD/m neutral density filters. A 5 kW methane burner flame was used 
to obtain the calibration factor. The heater was positioned such that the exposed surface 
of the material under test was 25 mm away. Finally, the system was checked for 


consistency by burning a 25 mm thick sample of black PMMA. 


Each time the incident heat flux was changed, the output from the heater was verified 


with one of the calibrated total heat flux gauges. 


1.4 Burning rate and ignition tests 

1.4.1 General 

The wood samples were provided such that the grain was parallel to the incident heat flux 
(i.e. cut across the grain) and perpendicular to the incident heat flux (i.e. cut along the 


grain) as shown in Figure 4. 
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Incident heat flux Incident heat flux 
parallel to grain perpendicular to grain 


50 mm 


96 mm 


Across grain (X) Along grain (L) 


Figure 4. Sample grain configurations. 


Four species of wood were tested in the study: Douglas fir, Redwood, Red oak and 
Maple. Douglas fir and Redwood are both softwoods whereas Red oak and Maple are 
both hardwoods. The samples were all cut from the sapwood portion of sections of 
lumber. The test numbering system used in this study was as follows 

nssot [i] e.g. IDFL2 
where 

n denotes the test series; 1 or 2 


SS denotes the species; Douglas fir (DF), Redwood (R), Red oak (O), 


Maple (M), 
O denotes the grain orientation; L is along or X is across (Figure 4), 
t is the test number for that species and orientation; 1 ... 7, 


and the optional post-fix letter ‘1’ indicates an ‘ignition only’ test. 
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Samples were stored in a desiccator at nominally 50 % relative humidity and 20 °C. The 
moisture content of each sample was measured with a hand-held moisture meter prior to 
exposure. In many cases, it was found that the moisture content of the samples was below 
the sensitivity of the meter. The minimum sensitivity value was recorded in such 


instances. 


All samples were tested in the horizontal orientation. Samples were wrapped in a single 
layer of aluminium foil, placed into the sample holder and backed by non-combustible 
ceramic fibre insulation material. In most tests the layer of ceramic fibre blanket was 
necessarily thin since the maximum height of the sample retainer frame is 50 mm and the 


nominal height of the samples was also 50 mm. 


The doors to the chamber were closed during the experiments and air was provided by a 
vent in the base of the chamber below the load cell. The spark igniter was located 
~10 mm above the surface of the sample. Sustained ignition is defined as when the 


sample continues to flame for an uninterrupted period of at least 10 s. 


The majority of the experimental tests used in this study were conducted at the University 
of Maryland by the author on behalf of Schroeder [10] as part of his analysis of the 
change in the structure of materials when exposed to an external heat flux for relatively 
prolonged durations. The test materials and experimental protocol were thus primarily 
selected by the requirements of Schroeder’s analysis. As it will be described in § 4.2, 


these requirements presented certain problems and limitations on the results obtained. 
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The experimental series’ were supplemented by a small number of additional tests that 
were not part of Schroeder’s study. The small number of extra samples available to the 
author limited these tests. The experiments described in this study were conducted in two 


series as described below. 


1.4.2 Burning rate tests 

The main ‘burning rate’ series of 54 tests included the complete measurement of time to 
ignition, mass loss, rate of heat release and smoke extinction data. Incident heat fluxes of 
25 kW/m’, 35 kW/m’, 50 kW/m’ and 75 kW/m’ were selected for these experiments. 
These heat fluxes are typical of values used by other researchers. For the majority of the 
burning rate tests, exposure times ¢, of 25 minutes were used, however, in a few cases the 


exposure time was extended to 75 minutes. A summary of the exposure conditions for 


each species of wood in the two grain orientations is given in Table 2. 


3 (2L, 1X) 


2 (IL, 1X) 
5 (3L, 2X) 6 (3L, 3X) 11 (6L, 5X) 


71, 4X) 6 GL, 3X) oe 25 (I2L, 13%) 


1(1L) 


Table 2. Main ‘burning rate’ series test exposure conditions for each species of wood in 
the two grain orientations. 


Each sample was nominally 50 mm thick and 96 mm square (measured to the nearest 


millimetre). For the majority of the “burning rate’ tests, 0.813 mm (0.032") diameter 
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sheathed Type K thermocouples were inserted into the samples through holes at heights 
of 4mm, 12 mm, 24 mm and 36 mm below the top surface (Figure 5). A thermocouple 
was located between the back of the sample and the retainer frame and a thermocouple 
was used to record the ambient temperature outside the Cone Calorimeter. The 
thermocouple measurements were recorded by a data logging system separate from the 


Cone Calorimeter data acquisition system. 


Cone heater 


Wood sample Spark 
a igniter 

° 

° 

fo) 

fe) 


Ceramic fibre 


blanket backing Retainer 
Thermocouples at 4 mm, frame 
12 mm, 24 mm, 36 mm 
and at back face 


Figure 5. Cone Calorimeter sample configuration in main ‘burning rate’ tests. 


Although the post-processing program provided with the Cone Calorimeter provided an 
average heat of combustion, in this study the value was obtained manually from the 
instantaneous data. This was necessary for several reasons. Firstly, at the end of each test 
the Cone Calorimeter normally requires that the sample remain on the load cell for 
approximately one minute while the data collected is completed. However, since the 


requirement of Schroeder’s work was to have exact exposure times, the sample was 
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necessarily removed immediately after the end of the test was declared. This resulted in 
the output from the load cell being affected and thus the final automatic average heat of 
combustion calculation. Secondly, the wood samples were exposed for long periods of 
time during which burning did not occur. As will be described in § 4.7, during these long 
test periods the instantaneous heat of combustion would shift as the wood charred and 
this char oxidised. Finally, as will be discussed later in § 4.2, a meaningful measurement 


of the instantaneous heat of combustion was not always achieved. 


Thus, where applicable, the average heat of combustion was found by obtaining the mean 
of the instantaneous heat of combustion data over the portion of the test between ignition 
and before any fluctuations at the end of the test (Figure 6). Where char oxidation effects 


were observed, a separate average heat of combustion for the char was obtained. 
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Figure 6. Obtaining the average heat of combustion (Test 1RL3). 


Observations were made prior to and during the burning of the sample. It was noted that 
as the samples burned, the flames would initially be at their tallest (Figure 7 {a}), then 
reduce in size over time and would sometimes disappear completely leaving the char 
layer to glow (Figure 7 {b}). Figure 8 shows an example of a normalised flame height for 


Test 2DFX9i (Douglas fir at 35 kW/m’) 
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Figure 7. Sample burning in the Cone Calorimeter; {a} 1:40 after ignition; {b} 25:00 
after ignition (Test 2DFX9i). 


0.8 


0.6 


0.4 


Normalised flame height [-] 


0.2 


0.0 
00:00 01:00 02:00 03:00 04:00 05:00 06:00 


Time [mm:ss} 


Figure 8. Normalised flame height (Test 2DFX9i). 


The flames were more likely to cease at the lower incident heat fluxes but that was not 
always the case. In some cases, after the majority of the exposed surface had ceased 
flaming, some smaller flames could be seen around the edges of the sample. This was due 


to shrinkage and volatile gases escaping from the vertical edges of the sample. 


1.4.3 Ignition tests 

A total of 41 additional ‘ignition only’ tests at lower heat fluxes were conducted where 
only time to ignition was measured. These tests were conducted between heat fluxes 
below 25 kW/m? down to heat fluxes in the region of the critical heat flux for ignition for 
a particular species of wood at a particular grain orientation. The critical heat flux is 
defined as the minimum external heat flux required to achieve piloted ignition of an 
exposed sample. In these ‘ignition only’ tests, the sample was exposed to the external 
heat flux until sustained ignition occurred or until it was determined by observation that 
ignition was unlikely to take place. If ignition occurred, the sample was extinguished 


immediately. 
The selection of a ‘failure to ignite’ criterion is somewhat subject to operator 
interpretation and patience. The ASTM standard for the Cone Calorimeter [7] suggests in 


paragraph 11.2.8 


If the specimen does not ignite in 10min, remove and discard, unless the specimen 


is showing signs of heat evolution. 


Pa 


Similarly, the ASTM standard for the Lateral Ignition and Flame Test (LIFT) apparatus 


[11] suggests in its paragraph 11.2.8 


The test is considered complete if ignition does not occur within 20 min. However, 


this is an arbitrary cut-off, and longer times can be considered. 


Clearly both of these test methods leave the ultimate decision as to when ignition has not 
(or will not occur) to the operator and the requirements of the particular experiment. In 
this study the decision as to when to terminate a test was of particular importance in 


determining the critical heat flux for ignition (§ 3.2.2). 


In the ‘ignition only’ tests single 50 mm thick samples of wood were cut into four equal 
thickness slices. For each test, the four slices were stacked in the sample holder to mimic 
the full thickness samples used in the main test series. The arrangement of samples in the 
Cone Calorimeter was thus similar to Figure 5 with the exclusion of the thermocouples. 
Critical heat flux measurements were not conducted for Red oak since all samples were 


used in the main ‘burning rate’ experiments. 


In retrospect the use of the stack of samples may have been unnecessary and also may 
have influenced the results somewhat. Where the stacked samples were subject to 
incident heat flux for a considerable period of time, the moisture and volatile content in 
the lower samples may have been altered so that when those samples were repositioned at 


the exposure surface in subsequent tests they were no longer effectively ‘virgin’ wood. 
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On the other hand, since samples were exposed for several tens of minutes before 
ignition, the moisture and volatile content would have undergone charges regardless of 


any effects. 


Prior to ignition it was noted that some samples would warp either away from or towards 
the cone heater. The warping was seen to ‘self-correct’ (i.e. return to almost level) in 
some instances. In addition, coupled with shrinkage, there were cases of the sample 
warping out of the retainer frame at one corner or along an edge. These factors may have 
introduced some margin of error into the ignition results since the sample may have 


ignited sooner or later than if it had remained level. 


In the tests conducted in this study it was found that the wood continued to ignite even at 
very low incident heat fluxes i.e. below 10 kW/m’, which is considerably lower than 
values quoted in the literature (see § 3.2.2). It was observed that at these low heat fluxes, 
a localised glowing could be seen on the surface of the wood prior to ignition. In such 
cases, ignition would eventually occur with the flames initially limited to the region of 
glowing but gradually ene over the exposed surface of the sample. In contrast, at 
higher heat fluxes the sample would immediately ignite over the complete surface of the 
sample. It is likely that this localised glowing contributed an additional source of energy 


to that provided by the heater to the surface of the wood eventually leading to ignition. 
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For example, Figure 9 shows the glowing on the surface of a sample of along grain 
Maple exposed to an incident heat flux of 12 kW/m’ (Test 2ML41). The glowing 


appeared approximately 1 hr into the test followed by ignition at 1 hr and 10 minutes. 


Figure 9. Localised glowing region observed on sample 2ML4i. 


Martin [12] alludes to this change in the ignition mechanism at low heat fluxes. He 
suggests that the ignition behaviour of cellulose can be split into three regions; 
convection-controlled, diffusion-controlled and ablation-controlled. He further notes that 
cellulose exhibits basically two kinds of gation phenomena without the presence of a 


pilot flame - spontaneous flaming and glowing ignition. 


The glowing ignition typically occurs under conditions of slow heating and therefore is 


not the usual response to the diffusion-controlled region. Its occurrence is dependent on 


pe 


the mineral content of the specimen and size of the irradiated area. Finally, the glowing 
ignition, almost without exception, sustains itself once started whereas the flaming 


ignition may be either transient or persistent. 


1.5 Literature review 

1.5.1 General 

There is a substantial volume of work in the literature regarding the pyrolysis, burning 
and charring behaviour of wood (and cellulosic materials). It is not within the scope of 
this work to be able to reference and review every study conducted. The exclusion of a 
particular piece of work does not neglect its importance to this field of study but simply 
highlights the inability of the author to cover every contribution. However, several of the 
more directly relevant studies have been selected as sources of material property data, 


modelling techniques and for comparison purposes. 


Kanury [13] gives a general overview of the ignition of solids by thermal radiation or 
convection. Roberts [14] reviewed the role of kinetics for the pyrolysis of wood and 
related materials. Simms [15] examined the role of thermal radiation on the damage to 
cellulosic solids by considering the chemical and thermal histories of the material. Work 
on char rate in wood includes studies by Kanury [16] who examined the phenomenon 
using Arrhenius pyrolysis kinetics. Schaffer [17] also investigated the charring rates of 


various species of wood. 
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Atreya and co-workers have done extensive work on the ignition and burning of wood. 
Included in this was the initial work by Atreya [18] that included experimental 
observations for the piloted ignition of wood and identified several important factors. A 
detailed finite difference ignition model was developed by Tzeng & Atreya [19] in an 
attempt to quantify the experimental observations. Abu-Zaid & Atreya [20] considered 
the effect of moisture on the ignition of cellulosic materials. Further work by Atreya, 
Carpentier & Harkleroad [21] examined the effect of sample orientation on piloted 


ignition and flame spread on wood. 


A detailed study of the pyrolysis kinetics of cellulose has been conducted by Suuberg, 
Milosavljevic & Lilly [22]. Cellulose was chosen because it is relatively well 
characterised material and, as described in § 1.2.1, cellulose forms a major component in 


the structure of wood. 


1.5.2 Ignition and burning rate models 

Several models for the burning rate of solid materials, both charring and non-charring, 
have on developed. Examples include the studies by Delichatsios & de Ris [23]; Chen 
et al. [24]; Wichman & Atreya [25]; Yuen et al. [26] and Parker [4]. These models range 
from simple treatments of the ignition and burning process, integral model approaches 
through to an analysis of the complex chemical kinetics involved in the burning of a 
charring material. Many of the models consist of complex computational codes that 
require a relatively large number of property values to complete their predictions. These 


factors can (at least at present) limit the use of such models since many of the properties 


yay) 


are difficult to obtain and the codes may not be suitable for incorporation into more 


general fire hazard models. 


In this study, we examine the integral model initially developed by Quintiere [27]. A one- 
dimensional pyrolysis model which includes the processes of charring, vaporisation, 
flame and heat conduction effects was proposed. This model was further developed by 
Quintiere & Iqbal [28] to solve the one-dimensional unsteady heat transfer equations 
during the pre-heating and gasification periods using an integral method. Anderson [29] 
studied the integral solution to Quintiere's [27] model and compared the integral solution 
with the exact solution. Finally, in the study conducted by Hopkins [30], the model 
developed by Quintere [27] was compared against experimental data for non-charring 


thermoplastics tested in the Cone Calorimeter. 


The integral model approach to the burning of a charring material has been successfully 
demonstrated by Moghtaderi, Novozhilov, Fletcher & Kent [31] using a computational 


fluid dynamics model. Simulation times of about eight hours are quoted by the authors. 


1.5.3 Experimental data 
Experimental data for the ignition and burning rate of several species of wood have been 


presented in the literature. Studies of particular reference are outlined below. 


Janssens [1], [32], [33] tested several species of wood in the Cone Calorimeter of which 


Douglas fir and Redwood will be of particular interest to this study. In all cases the 
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samples were tested in the vertical orientation with the grain perpendicular to the incident 
heat flux (i.e. equivalent to the along grain orientation used in this study). Parameters 


obtained included critical heat fluxes, ignition temperatures and energy release rates. 


Measurements of the burning rate of wood have been made by Tran & White [34] using 
the Ohio State University (OSU) apparatus. Four types of wood, namely Redwood, 
Southern pine, Red oak and Basswood were tested at a range of incident heat fluxes 
between around 17 and 56 kW/m’. Parameters including times to ignition, average 
burning rates, charring rates were obtained. The study also obtained properties for the 


char residue for each type of material tested. 


Other experimental data have been reported by Parker [35] for Douglas fir and 
Dietenberger [36] for Redwood in both the LIFT and Cone Calorimeter. The study by 
Hopkins [30] also includes data for two charring materials, namely Redwood and Red 
Oak, but no detailed analysis or comparison with the integral model was conducted with 


these data. 
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2. THERMO-PHYSICAL MATERIAL PROPERTIES 

The ignition and burning rate integral model used in this study requires a number of 
properties to be obtained for the material. A few of the properties can be easily measured, 
others can be obtained from experimental data and the remainder may be obtained from 


the literature. 


The determination of fundamental ciate properties can be a complex process. For the 
pyrolysis of wood we require the thermal conductivity k, density p and specific heat 
capacity c and the related properties of thermal inertia kpc and thermal diffusivity k/pc. 
These properties may change as the material undergoes thermal and/or mechanical 


changes. 


In addition we may also need to obtain estimates of the heat of gasification of the 
pyrolysis gases L and heat flux of the flame qj, as the material burns. Anderson & 
Quintiere [37] suggest that an iterative procedure of choosing a heat of gasification and 
comparing the predictions of the model with experimental data be used. The guess is 
refined until there is good agreement. This approach is used in this study as the way in 


which to obtain reasonable values for the heat of gasification and flame heat flux. 
2.1 Density 


The density of wood is primarily dependent on the species but it will also vary by 


individual tree and within that individual tree. Any moisture in the wood will also affect 
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the density. In this study, the average bulk density of each sample was calculated from its 


mass and volume such that 


M,, 0 


Pee oom 
A, Myo 


Equation 2 
where Ap is the surface area of the exposed face of the sample and 7,,, 9 is the initial height 


of the sample. In this study, the mass and volume were recorded prior to testing in the 
Cone Calorimeter and the density simply obtained from Equation 2. The calculated 


density for each individual sample is shown in Figure 10. 
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Figure 10. Sample density variation. 


The sample densities for Douglas fir exhibit a spread with the along and across grain 
samples clustered into two distinct groups. Similarly for the Red oak samples. The 
Redwood samples also show a spread of values but no apparent clustering. The densities 
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of the Maple samples show less variation than found for the other three specie. Table 3 
shows the average density of the wood species tested with a comparison of some values 


quoted in the literature. 


4301] 


660 [34] 


465 [1], 500 [35], 450 [2], 
512 [40] 


Table 3. Average density of virgin wood samples. 


2.2 Specific gravity 

The specific gravity s of wood is based on its weight when oven dry and volume at 6% 
moisture content [38]. The relationship between density and specific gravity, using values 
for the specific gravity quoted in the Wood Engineering Handbook [38] for the species 


tested in this study and the average densities given in Table 3, is shown in Figure 11. 


af 


0.70 


0.65 . Maple % 
0.60 esc ve Se We a A 
= 055 note 
sare Pat 
2 i oe 
@ 0.50 Rick 
tT) nee 
a - 
S aS 
2 0.45 og ~~ Doughas fir 
i? a] are 
0.40 afte 
ae Closed symbols - along grain 
@ _-O Redwood Open symbols - across grain 
0.35 


300 350 400 450 $00 $50 600 650 700 750 800 
Average sample density, p, [kg/m] 


Figure 11. Relationship between density and specific gravity. 


It is clear from Figure 11 that the measured average densities of the samples are 


consistent with the average specific gravity values. 


2.3 Moisture content 

The moisture content of wood may be assumed to be a pseudo-property of the material 
and, as described in § 2.4 and § 2.5, it can have an influence on the thermal conductivity 
and specific heat capacity and thus the ignition and burning rate characteristics of wood. 
The moisture content is a function of the species of wood and the conditions in which it is 
stored (§ 1.2.1). For example, Cholin [3] shows that the moisture content of wood is 
affected by the relative humidity of the atmosphere. The study by the Fire Officers 
Committee quoted by Cholin [3] demonstrates how the pers in the moisture content 


of wood increases the time to ignition for a given incident heat flux. 
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The measured moisture content of the samples tested is shown in Figure 12. The 
minimum sensitivity limit of the moisture meter (see § 1.4.1) for each species is also 
shown. Thus where the moisture meter indicated that the moisture content was below the 
minimum sensitivity it was no possible to determine the exact moisture content but only 


determine the maximum limit. 
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Figure 12. Measured moisture content of wood samples used in "burning rate' tests. 


For the Red oak and Maple samples, the moisture content was at or just above the 
sensitivity limit of the meter. The majority of the Douglas fir and Redwood samples had 
moisture contents above the sensitivity limit of the meter with the Redwood exhibiting 


the widest variation in moisture content. 
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2.4 Thermal conductivity 

The study by Fredlund [39] showed how the thermal conductivity varies in wood with 
emittance, density, moisture content, temperature and the type of gas enclosed in the 
material. Thermal conductivity increases markedly with increasing moisture content, 


being about twice as high at 100 per cent moisture content as it is at 10 per cent. 


The thermal conductivity also depends on the orientation of the grain of the wood. 
According to the Wood Engineering Handbook [38], the thermal conductivity of wood is 
approximately 2.0 to 2.8 greater along the grain than perpendicular to the grain. Fredlund 


[39] quotes a study that gives the range of ratios as between 1.75 and 2.25 (though he 
appears to use a weighting factor of ye 537 1.72 in his calculations). Desch & 


Dinwoodie [40] quote values for the thermal conductivity of Spruce and European oak 
for the various grain orientations (parallel or tangential/radial). The ratios of the values 
give 2.10 and 1.75 for Spruce and European oak respectively. From these literature data it 
is assumed that the increase in the thermal conductivity for the samples tested across the 


grain is typically 2.1 times greater than along the grain for any species of wood. 


Janssens [41] provides equations that allow the thermal conductivity to be calculated. The 
equations include the effect of moisture, acne the thermal conductivity contributions 
due to the air and radiation effects. Figure 13 shows typically how the thermal 
conductivity for dry and wet wood varies with temperature and density using Janssens 


equations. 
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Figure 13. Variation of thermal conductivity with temperature for typical wood densities. 


Many investigators have suggested that the thermal conductivity of dry wood varies as a 
linear relationship with temperature and Figure 13 shows this relationship. Figure 13 
indicates that the inclusion of moisture results in a non-linear relationship although it is 
not clear over what range of temperatures Janssens equations are valid. Finally, using 
Janssens equations, the thermal conductivity for wood with a density of 500 kg/m’ at 

50 °C was calculated to be 0.13 W/m.K with 10% moisture and 0.36 W/m.K with 100% 


moisture. This result is comparable with the increases suggested by Fredlund [39]. 


2.5 Specific heat capacity 
The specific heat capacity of wood depends on temperature but is practically independent 


of density or species. For oven-dry wood, Desch & Dinwoodie [40] give the specific heat 
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capacity as 1,360 J/kg.K. When wood contains water, the specific heat is greater than dry 
wood because of the larger specific heat of water. The apparent specific heat of moist 
wood is larger than the simple sum of the separate effects of wood and the water. This is 
due to the thermal energy absorbed by the wood-water bonds. The specific capacity of 


wood will also increase with increasing temperature. 


2.6 Char fraction 

As described in § 1.2.2, the burning of a charring material consists of a pyrolysis front 
moving through the unburned material leaving behind the char residue. Figure 14 shows a 
schematic of the form of a wood sample prior to testing and after exposure to the incident 
heat flux. Down to a certain depth below the exposed surface of the sample becomes char 
and the top surface of this char layer 71s lower than the original sample height 779. In 
practice, the top surface of the char layer and the interface between the char layer and the 
virgin wood below were not necessarily horizontal and the interface not necessarily as 
distinct as suggested in the schematic (see § 4.5.2). Thus the determination of the final 


height of the sample and the char depth were open to a certain degree of interpretation. 


At the termination of a test in the Cone Calorimeter, the remainder of the sample was 
weighed and its dimensions recorded. An average final height and char layer depth were 
reported by Schroeder [10] and those values are used in this study to obtain the char 


fraction as follows: 
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Figure 14. Pre-test and post-test char fraction determination. 


Post-test, the mass of the whole sample is 
thus the mass of the char is 


LPS Mal 9 sei Pay 


The mass and volume of the post-test virgin wood can be found from 
my ¢ = Pw Vos 


yay re A,(7, —15,7) 


thus 


Mm, ¢ = Pw,s Ay oF 1,7) 


Substituting into Equation 3 gives 
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Equation 3 


m,=M,— Py f A,(n, —15,7) 


hence 


Vo. Pp = My — Pug A(T — 1g.) 


Equation 4 
The volume of the char portion is 
Vor = Ants 
thus the density of the char, from Equation 4 
_ Ms — Pug Ay —M,7) 
Pi atamnaaan cen aaa 
Ay 1g 4 
The char fraction is defined as 
Ppuesd £4 
gat 
Pw 
thus 
ie aie mM, — Py A,(1; a 1,7) 
Pw A, Ny f 
Equation 5 


and it is assumed that p,, , = /,,o, i.e. the density of the post-test virgin wood is the same 


is the initial density of the sample. 


At temperatures exceeding 200 °C, the thermal degradation process results in a decrease . 
in the volume and the amount of shrinkage is different parallel or perpendicular to the 
grain. Janssens [41] provides methods to obtain the shrinkage as a function of 


temperature and wood type. Analysis of the post-test samples by Schroeder [10] 
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identified varying amounts of shrinkage compared to the original sample. The shrinkage 
factor ® was used to obtain a modified char fraction such that 


$= af 


Equation 6 
Tran & White [34] suggest a similar relationship between the char contraction factor and 


the char yield in their analysis. Thus the char fraction can be obtained from 


a a —- Pp, A(n, aod 


AN, + 


Equation 7 
The char fractions obtained from the test samples are detailed in § 4.5. 
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Figure 15. Thermal conductivity of dry wood and char with density. 


Figure 15 shows the variation in the thermal conductivity of a typical dry wood at given 


densities as given by Janssens [41] (§ 2.4) and the variation in the thermal conductivity of 
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char [41]. Figure 15 shows that the change in thermal conductivities is almost equivalent 
over a wide range of densities. As discussed in § 2.1, the average density of the samples 
tested in this study ranged from around 330 kg/m’ up to around 750 kg/m°. Thus, it can 
be assumed that the ratio of the thermal conductivity of the virgin wood with the thermal 


conductivity of the char can be approximated to the ratio of the virgin wood density to the 


char density such that 
Ky Po 
soi of te 


Equation 8 
Figure 16 confirms that Equation 8 holds for a range of virgin wood densities with the 


approximation better at the higher virgin wood density of 800 kg/m compared to 


300 kg/m’. 
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Figure 16. Thermal conductivity and density ratio approximation. 
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2.7 Thermal inertia and thermal diffusivity 
The thermal inertia J is defined as the product of the thermal conductivity, density and 
specific heat capacity such that 

I=kpc 


Equation 9 
and the thermal diffusivity @ is defined as 


Equation 10 
The integral model described 1n this study requires values for the thermal conductivity k 


and the specific heat capacity c at elevated temperatures. It is clear from § 2.4 and § 2.5 
that, since the thermal pear and the specific heat capacity are temperature dependent, 
the thermal inertia at ignition is not that what would be obtained at ambient conditions. 
Instead the thermal inertia at ignition is an ‘apparent’ value and it will be shown in § 


3.2.4 that this apparent thermal inertia can be obtained from ignition data. 


In the study by Parker [35] it was shown that the thermal diffusivity of Douglas fir 
remained at an almost constant value of 2.1 x 10°’ m’/s up to temperatures of around 
250 °C. Similarly, Suuberg et al. [22] found that the thermal diffusivity of raw cellulose 
remained constant at 0.86 x 10°’ m?/s +22% between 116 °C and 289 °C. Janssens [41] 
also quotes work in which it is suggested that the thermal diffusivity is independent of 
temperature. Thus, in this study it is assumed that the thermal diffusivity is constant for 


each given species of wood. 


42 


The values of the thermal diffusivities perpendicular to the grain used in this study are 
given in Table 4. The value for Douglas fir is taken from Parker [35] and the value for 


Redwood taken from Dietenberger [36]. 


Table 4. Thermal diffusivity values. 


Values for Maple and Red oak were not found in the literature and so these had to be 
estimated by simply taking the average of the two known values quoted. Taking an 
average value was justified by the fact that according to the Wood Engineering 
Handbook [38] the typical value for the thermal diffusivity of wood is 1.61 x 107 m’/s 
and this value varies with specific gravity s over the range of 0.35 to 0.60 by 

6.45 x 10° m’/s. Assuming that the typical quoted value in the Wood Engineering 
Handbook for the thermal diffusivity is at the mid-range of the specific gravity (i.e. for a 
specific gravity of 0.48), the variation of thermal diffusivity with specific gravity can be 
compared with the values quoted by Parker and Dietenberger assuming similar moisture 


contents for each species (Figure 17). 
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Figure 17. Variation of thermal diffusivity with specific gravity. 


It is clear that the data quoted by Parker and Dientenberger closely match the relationship 
given in the Wood Engineering Handbook. Furthermore, there is almost negligible 
change in the thermal diffusivity over the specified range of specific gravities. Thus, 
given the relationship between the density and specific gravity discussed in § 2.2 (Figure 
11), using an average value for the thermal diffusivity for the Maple and Red oak appears 


reasonable. 


It is interesting to note that the thermal diffusivity of 6.6 x 10°’ m’/s used for White pine 


by Moghtaderi et a/. [31] was considerably greater than the values quoted above. 


Using Equation 9 and Equation 10, apparent values for the thermal conductivity and the 


specific heat capacity can be obtained. However, as discussed in § 2.4, the thermal 
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conductivity of wood increases by a factor of around 2.1 along the grain compared to 
perpendicular to the grain. Hence a factor fof 2.1 is introduced for the across grain 


orientation thermal diffusivity. Thus, from Equation 9 


pe= 


TN 


Equation 1] 
substituting into Equation 10 and including f where f= 1 for the along grain samples and 


f=2.1 for the across grain orientation, 


Equation 12 


rearranging to solve for k we get 


k? 
a=— 
fl 
Equation 13 
k= Jfal 


Equation 14 


and we can find c from 


Equation 15. 


2.8 Emissivity 
In this study the emissivity of the surface of the wood is assumed to be unity. Janssens [1] 
quotes several sources regarding the emissivity and absorptivity of wood that suggest the 


absorptivity of virgin wood is on average 0.76, independent of species. After thermal 
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exposure begins, this value changes due to the darkening of the surface as it chars. The 
emissivity of oven dry wood varies between 0.60 and 0.72 depending on species. Finally, 
the assumption that Kirchoff’s law (a@ = €) holds is reasonable for most analyses. Thus we 
might expect the emissivity to be around 0.72 prior to exposure and this value to 


approach 1 as the surface chars due to the external heat flux. 


46 


3. IGNITION MODELLING AND EXPERIMENTS 


3.1 Theory 
3.1.1 The integral model 
The integral analysis for ignition was developed by Quintiere [27] assuming radiative 


ignition of the surface of the material due to an incident radiative heat flux. 


bE? 
bee: 


Figure 18. Integral model ignition scenario. 


The following assumptions are made for the ignition model: 
(a) Ignition occurs when the surface temperature achieves a critical value, Tig — 
(b) Solid is inert up to ignition and 


(c) Solid is infinitely thick 
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The thermal heating of the solid is depicted by a thermal penetration layer of depth d(t) as 
shown in Figure 18. By considering the incident heat flux and the losses due to radiation 


and convection, the net heat flux at a given time ¢ is given by 
q"(t)=adj-eo(T; -T')-h (I, -T) 


where the emissivity ¢ and the absorptivity a are assumed to be 1 as discussed in § 2.8, 


thus g"(t) becomes 


q'(t)=4[-o(T! -T;) -h,(L, -T) 


Equation 16 


The control volume form of conservation of energy 
— alle + Jone - w).ndS = ORS: nal dAeior ys [pw. ndS 


Equation 17 


gives 


pox T(x)dx+ pcT, Gar q'() 


p| = [ 7Q)ax-r, a). q’(t) 


Equation 18 


48 


Where c is the specific capacity of the wood at ignition i.e. c=c,, ,, and pis the density of 
wood at ignition i.e. 9 = P,, ;, = P,,Since we assume the density remains constant. 


Consider 


< f(@-n)de 


Equation 19 


=< [rax-< [Tae 


_d fs iba ps 
wait , Tax-T— |e] 


ea [Tax-7, 2 
dt *° dt 


Thus, from Equation 18 and Equation 19 
d es 
pe [(T-T)ar= 4") 


Equation 20 


We need to select a temperature profile through the region 5 such that the boundary 


conditions are 


oT 
henx=0, q’(t)=-k— 
(i) when x q (t) An 
(i) whenx=6, T=T, 
(iii) when x=6, ae = i.e.no heat loss 


where k is the thermal conductivity of the wood at ignition i.e. kK =k,, ,, . Consider a 


profile such that 
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Equation 2] 


when x = 6 


‘which satisfies boundary condition (11); by expansion of Equation 21 


on 2 
r-7,=4 (20), _ 2% x 
2k 


and differentiating both sides with respect to x 


<(r-1,)= BUD of. ae + x ) 


DOP Or eieroo On 


or he = £08 2+) 


OF Ane 
x  k  ké 


when x = 0 


i A) 
k 
oT 


mhesy Ae) 


which satisfies boundary condition (i); when x = éd then 
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ar _ FO FOS 
Ox k ké 


which satisfies boundary condition (111). Substituting Equation 21 into Equation 20 gives 


d pq°(t)d 
dt” 2k 


dt” 2k 


{i- x) a= q'(t) 


gO 5 s{- 2 x 
5 8 


= Jar q'(t) 


on 3 5 
pi £0! a -x24 | =4") 


dt 2k 


dime Kis | widen ae 


pb £01 55 oO + 


2 FO 


=@q'"(t 
PAGES Ae 


d k 
—q"(t)d? =6—q"(t 
ow: pao 


Integrating 


at 
= 6— | g"(t)dt 
bt 


Zh 70) 


ant 


Figure 19. Assumed heat flux profile. 


From Figure 19, we assume that the incident heat flux is the average of the heat flux at 


time = 0, i.e. g”(0) and at time =t, i.e. g’(t) 


[a@at = pe ££) | 
thus 

PLAC KOET KO 

gr"(t)d -6£| ; | 


Equation 22 


From Equation 16, at time ¢ = 0, the surface temperature is at ambient i.e. JT, = 7, thus 
q'(0)=47 


Equation 23 


and at time ¢ = tig, we assume that the surface temperature is at the ignition temperature, 


T, = T,, thus 


ae 


q(t, )=4f-—o(T! - Te) -h,(T, -T,) 


Equation 24 
From Equation 16, multiplying through by -1 


-q'(t)=-q?+o(T! -T) +h (I, -T) 


and then dividing through by —q7 . 


Gt) | SG sith. Gat) 
q q 


Let 


o(f' ~ Ti) +h(T,-T,) 


B= - 
qi 


Equation 25 


FO) 
qi 


q’(t)=4701- B) 


Equation 26 
The parameter £ characterises the magnitude of radiation and convective losses relative 


to the incident heat flux. Substituting Equation 23 and Equation 26 into Equation 22 we 


obtain 


Pays ii Aan 
qi(l f)6 eal > | 


past Made 
pe| qi(1- B) 


53 


paak|28 
pc, 1-B 


Equation 27 
If we consider the surface at the time of ignition and assume that the surface temperature 


is at the ignition temperature t = fig, YT, =T, , x = 0, then from Equation 21 


Equation 28 
where from Equation 26 
q't,,)= ql ie B.) 
Equation 29 
and from Equation 25 
ga 2a ~ Te) + hy ~B) 
i qi 
Equation 30 


When 9; is large, from Equation 25, ig > 9, thus from Equation 28 
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oy, peu ae 


3" 2-0 |[Gt—OF 


Equation 3] 


The 2/3 coefficient has been found to be 7/4 in the more exact (pure convective loss) 


solution of this problem (see § 3.1.2). Substituting Equation 29 into Equation 28 we 


obtain 
4 1- B, ig (r,, ar Tal 
lig kpc ond 
3 2-f, qi la De 
4 (7, -7,F 
lig oa Eee 2 1 Py) 
am Biz aye B. qj 
Equation 32 
Or alternatively 
(r, - 7) 
lie *s C,,kpc = ; wn 
Equation 33 
where 
5 2-2, 0-8) 
Equation 34 


As £,, 1, Equation 34 approaches oo and thus from Equation 33, t,, 0 also. As the 
time to ignition increases we are approaching the critical heat flux for ignition. From 
Equation 30 with £,, >1 


2B) 


_ oT -Ty) +h, Ty) 
qi 


] 


or alternatively, with Yo =9 as > © 


q:, = O(T, -T )+h(T, -h) 


Equation 35 
Thus from Equation 30 
24- 
Aer 
Equation 36 


3.1.2 Comparison of See atiiate solutions for ignition 

a) Quintiere [42] compared the approximate solution for ignition from the integral model 
with the exact solution for convective heat loss only and the approximate solutions by 
Delichatsios et al. [43]. For the exact solution for convective heat loss only, Drysdale [2] 


states that 


TT at _\_.{ Jat 
an tea ‘ ee 


Given, from Equation 16, for convection only 
q; = h, (7, oe T,) 
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thus 


T., -T. 2 
__° —=]-e” erfc(y) 
qi 
G 
where 
t. 
=} (az 
fldsleae 


From Carslaw & Jaeger [44], when y > 0 


Vx er 
— rf ee 
- e c(y) oy 


or rearranging Equation 39 


-y? 
erfe(y) = as 


Then, from Equation 37 and Equation 40 


rie: 
e’ 


lim(1—e” erfe(y)}=1-e” 


ae TES 


1 
=1-—— > 1 
Vary 


From Equation 35, neglecting radiation heat losses, as t,, > 0 


Equation 37 


Equation 38 


Equation 39 


Equation 40 


qi =h. Tig — To) = 4c, 


Equation 41 
When ¢,, > Oor vy + 0, from Carslaw & Jaeger [44], 
2 
erf (vy) = —= 
(y) RES Y 
Equation 42 
and 
erfe(y)=1—erf(y) 
Equation 43 
By expansion of the nght-hand side of Equation 37 
2 2) 2 
1-e” erfe(y)~1-(l+y +\1-er+..| 
Vx 
wz1—1+ mae y 
1 
ee 
Vx 
Equation 44 
Thus substituting Equation 44 and Equation 38 into Equation 37 
sig tao eee a eg 
G Un 'Vkpe 
h, 
Equation 45 


therefore, the time to ignition is 
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Equation 46 


To make the approximate integral solution given in Equation 32 fit the limit of the exact 
(pure convective loss) solution, let 


Z 


(2 a Biz tt aa Bi, ) 


he 
4 

Equation 47 
where Z is a new constant in place of the 4/3 given in Equation 32. At high incident heat 
fluxes, ¢,, >0 and f,, 0, thus 


== or Zee 
4 2 


NIN 


therefore, substituting for Z the 4/3 in Equation 32 we obtain 


= 
lin = koe ( . -1) 
2" 12-4,)I-4,)} 4 
Equation 48 


or 


1 1 


a —2.M-B, cs 
Ate It kpc 2\2 BN Biz T,, -T, 


b) In the study by Delichatsios, Panagiotou & Kiley [43], the authors suggest that when 
the incident heat flux is greater than about three times the critical heat flux (i.e. g7 > 3q7, ) 


then 


1 2 1 


= (-7) 7; — 0.644’ 
ti, lr kpc ip, ek 7. [4 q’,| 


Equation 49 


and when the incident heat flux is less than 1.1 times the critical heat flux (i.e. 


g; <1.19%,) then 


Equation 50 


where 7, is the pyrolysis temperature and the critical heat flux ignores convective heat 


fluxes which are considered by Delichatsios et al. to be negligible, that is 


Equation 51 
Defining the following dimensionless heat flux (from Equation 30) and time variables as 


on 2 
Lo eegn er fj ; 
—ett and 7, = (7, 1)’ kpc = respectively, we can write Equation 48 as 
B, ig Ver ‘ ip pv T, ; kpc 
1 


Equation 52 
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or by rearranging 


Equation 53 


plus the Delichatsios high and low heat flux equations (Equation 49 and Equation 50) as 


-7[7--064} Pyar, 


Equation 54 


Equation 55 


Figure 20 shows a comparison between the integral model and the Delichatsios high and 
low heat flux equations. The solutions to the two models run parallel at high heat fluxes 


and both models terminate at the same point at the intercept to the x-axis. The 
Delichatsios solutions are given for specified limits els (Equation 54 and Equation 55). 
ig 


By extending the two solutions for the region between the specified limits such that they 


overlap we find that they cross at around Beg 6. 
ig 
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Figure 20. Comparison of the integral model and the Delichatsios et al. equations. 


Figure 20 also includes the extrapolation of the high heat flux portion of the integral 
model which shows that there is an error in using such an extrapolation compared with 
the integral model solution for the determination of the intercept along the x-axis. Let 


TE areas be the intercept of the linear extrapolation of a graph of = plotted against 
. T ig 


B ig,intercept 


ie From the integral solution given in Equation 53, choosing values of sgh for two 
ig ig 


typical high heat flux cases we get 
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Equation 56 


By assuming that the difference between the exact solution and the extrapolated solution 


] l 
iS (4 7 a , we can write Equation 53 as 


i Die meron 


ees | alee sl sg, 
Tig Va By B, ig, intercept 


Equation 57 
Substituting in for our approximate values we obtain 


4.787 = (59 : a ee 9 758 
Va ig, intercept B, ig, intercept 
1.954 = F(2-g2 | Be 0768 
1 B ig, intercept B ig, intercept 
qi 


Hence the extrapolated intercept is g’ ~ 0.76 | thus the critical heat flux is found from 


cr 


Equation 58 
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In comparison, from Equation 53 (which is equivalent to Equation 57), Delichatsios et al. 
obtain a correction factor of 0.64 for the determination of the critical heat flux using high 


incident flux data. 


3.1.3 Analysis of the C;, parameter 

A qualitative comparison of the Cg parameter between the integral model and other 
solutions can be conducted. Equating the Delichatsios et al. [43] low heat flux equation 
(Equation 50) with the integral theory presented in this study (Equation 28) and assuming 
the pyrolysis temperature used by Delichatsios et al. is equivalent to the ignition 


temperature 7, = Tig we obtain 


Equation 59 
thus 


Equation 60 
Solving for fiz numerically we obtain a value of 0.67. Similarly, equating the 


Delichatsios et al. [43] high heat flux equation (Equation 49) with the integral model 


(Equation 28) we obtain 


(a eee ne rh) 640" | 
4 ta at Vakpe (T, ~T,) 
‘eal esre 
ig 


Equation 61 
thus 
Py eee lays 1 2 
*" SS OES SSS I q.) + 0.644—— 
mK pc ae l= Pe 4, lp. mk pc 
A wale 3 ert Wee Pi 
Equation 62 
At high heat fluxes g/ >> q?, thus —> — 0 and therefore 
rete piosccel ait Lynztsll ta 
 7kpc 1B. 
he Bie 
3 1 Big 
Equation 63 
or : 
Sie B. 32 
Zhe 16 
Equation 64 


However, this time solving for f,, numerically, we obtain a value of -0.433. Using the 


definition of Cjg given in Equation 34 such that 


4 1 
Cs e te a Big tt oat Bie 1 
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we obtain values for Cig of 3.24 and 0.38 for the low heat flux and high heat flux cases 


respectively. 


Abu-Zaid & Atreya [20] found experimentally that Ci, depends on g/ with values from 
0.5 to 0.62 obtained for heat fluxes of 10 to 60 kW/m’. Figure 21 shows the qualitative 
comparison between Abu-Zaid's results, the analysis of Delichatsios et al. [43] equations 
and the theoretical values for Cig using the integral model with a critical heat flux selected 


to be a typical value of 12 kW/m’. 


Best fit lines are shown through Abu-Zaid's measurements and the Delichatsios values to 
examine the trends in the relationships. At high incident heat fluxes, the integral model 
compares well with Abu-Zaid's data and Delichatsios equations. At low incident heat 
fluxes there is a significant difference between the Abu-Zaid's data value and the integral 
theory. The Delichatsios equations fall somewhere between the integral model md Abu- 


Zaid but with the trend at least similar to the integral model. 
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C, with 4/3 coefficient (1) 
C;, with 27/4 coefficient (2) 

Abu-Zaid & Atreya's experimental results (3) 
Inferred from Delichatious et al. analysis (4) 


Incident heat flux, q” , [kW/m’] 


Figure 21. Form of Cig parameter with incident heat flux. 


3.2 Analysis 

3.2.1 Times to ignition 

Time to ignition data for each wood specie in the across grain and along grain 
configurations are given in Table 5 to Table 8 and plots of ignition time against incident 
heat flux were made (Figure 22 to Figure 25). The plots include a ‘best-fit’ trend lines (in 
the form of a power curve, ax’ , where a and 5 are different constants for each set of data) 


through the data to indicate the form of the data. The data and ‘best-fit’ curves confirm 


that, from Equation 28, as qj — q?, then t,, «© and when q/ >> q7, then ¢,, > 0. 
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0.0110 12.0 
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0.0722 10 Be) 
0.0700 22.0 
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0.0020 
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0.0539 31.0 
0.0219 370.0 
0.0015 493.0 
0.0012 564.0 
0.0008 752.0 
0.0005 806.0 


0.0003 1685.0 | 0.0006 
0.0002 1465.0 | 0.0007 
2395.0 | 0.0004 

73.0 0.0137 


Table 5. Ignition data for Douglas fir. 
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Figure 22. Ignition time against incident heat flux for Douglas fir. 
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Figure 23. Ignition time against incident heat flux for Redwood. 
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Figure 24. Ignition time against incident heat flux for Red oak. 
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Table 8. Ignition data for Maple. 
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Figure 25. Ignition time against incident heat flux for Maple. 
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With the exception of the Red oak, it is noted that the across grain orientation samples 
ignited sooner at the lower incident heat fluxes (~20 kW/m’ or below). Yet at the high 
heat fluxes the along grain samples ignited quicker. Further analysis of this observation is 


given in § 3.4. 


3.2.2 Critical heat flux 

As described in § 1.4.3 the critical heat flux is taken to be the minimum incident heat flux 
required to obtain sustained piloted ignition of a material. The critical heat flux can be 
obtained experimentally by successively exposing samples of the material at decreasing 
incident heat fluxes until ignition no longer occurs. Thus the critical heat flux is 
somewhere between the lowest incident heat flux at which ignition occurred and the 
highest incident heat flux where ignition did not occur. Clearly this approach can be a 
time consuming process as it may require several tests to find the bounds of critical heat 
flux depending on the resolution required. In addition, as the critical heat flux is 
approached, then times to ignition become increasingly longer. Finally, there is the 
question as to how long one should wait before deciding that ignition will not occur. As it 
will be demonstrated in this study, ignition may not occur until anything between several 


tens of minutes and up to one and a half hours have elapsed. 


ths 


(*) no ignition case not achieved 


Table 9. Critical heat fluxes obtained from experiments. 


Table 9 shows the critical heat fluxes obtained aS the ignition experiments. It should be 
noted that in some cases the final experimental critical heat might not have been achieved 
simply due to the lack of test samples. In some cases, the sample at the lowest measured 
incident heat flux achieved ignition and thus the case in which no ignition was ever 
obtained was not found. Thus ignition may still have been possible at even lower incident 
fluxes. Such data are indicated in Table 9. The glowing ignition process at low heat 
fluxes was somewhat unexpected and had this been identified earlier in the study then the 
ignition testing protocol would have been modified accordingly in an attempt to gain 


further insights into the mechanism. 
As an alternative to directly obtaining the critical heat flux from an experimental 


procedure, the critical heat flux can be obtained from time to ignition data. By 


rearranging Equation 28, we obtain 


me 


tig Ae 1-B, | Vie -% 
2- B; 
Equation 65 
From Equation 24 and Equation 35 
ali, )= 41-4, 
thus 
l ie. 1 1 
== ¢-¢) ———— 
t iz 4 k 1- B., I, i 1) 
Wied hay i 
Equation 66 


The critical heat flux can be found from the plot of 1/ te against incident heat flux. By 


plotting a best-fit straight line through the data, the intercept yields a low estimate critical 


heat flux where, from Equation 66 


slope = 


4 
—k 
3 kee 


and 


intercept = —slope.q;, 


This follows because by the above theory the slope of the curve becomes vertical when 
1/ ie — 0 or £,, =1, while for higher heat fluxes the curve has a distinct positive finite 


slope as given by Equation 31 (see Figure 20). 
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In the study of non-charring materials by Hopkins [30], it was suggested that a linear 
regression through data below 40 kW/m’ gives a better measure for the critical heat flux 


since at lower heat fluxes ignition takes longer. However, examination of the plot of 
1/ (is against incident heat flux show that at low heat fluxes the data tends to exhibit a 


secondary upward trend towards a very low critical heat flux. This was as a result of the 


localised glowing ignition discussed in § 1.4.3. 


Thus, it was decided that the critical heat flux without the effect of the localised heating 
could be obtained from a linear regression through only the ‘high’ heat flux 
measurements and the low heat flux data was not utilised in the determination of the final 
critical heat flux. In this case, the ‘high’ heat flux data was taken to be where the incident 


heat flux was around 20 kW/m’ or above. The selection of the lower limit of the ‘high’ 
heat flux data was based on experimental observations, the shapes of the 1/ Af curves 


and from the theory shown in Figure 20. In this case ‘high’ heat fluxes are taken to be 


where q; > q/, and Figure 20 suggests that the integral model gives an approximately 


straight line when —— 21.5 i.e. gj 2>1.5q?,, in order to be consistent with the theory. 
ig 


Since values for the critical heat fluxes were found to be at most around 12 kW/m’ (and 


1.5 x 12 = 18) then a ‘high’ heat flux threshold of 20 kW/m?” is reasonable. 


Figure 26 to Figure 29 show the 1/ im against incident heat flux data for the four wood 


specie with a linear regression through the ‘high' heat flux points shown by large symbols. 
The figures also show the theoretical curves obtained from Equation 28 and the critical 
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heat fluxes obtained by the linear regression and theoretical curves. Table 10 shows the 


critical heat fluxes obtained from the intercept of the linear regression line. 
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Figure 26. Determination of the critical heat flux for ignition for Douglas fir. 
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Figure 27. Determination of the critical heat flux for ignition for Redwood. 
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Figure 28. Determination of the critical heat flux for ignition for Red oak. 
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Figure 29. Determination of the critical heat flux for ignition for Maple. 


In order to investigate the difference between using only the ‘high’ heat flux data and all 
of the data linear regression fits were also made through all of the time to ignition data 
obtained for the Douglas fir, Redwood and Maple species and the critical heat flux 
determined. These data are compared with the ‘high' heat flux data critical heat flux 
values in Table 10. In general, the along grain orientations show little difference. 
However, the across grain orientations for the Maple and Redwood show significant 
differences with the ‘high' incident flux data giving critical heat fluxes approximately 
twice are large. There is no difference between the data for the Red oak since no low 


incident heat flux measurements were made in the experiments. 
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Table 10. Comparison of critical heat fluxes for ignition using ‘high’ and all incident heat 
flux data. 


Table 11 shows the final critical heat fluxes for ignition for each specie using the 'high' 
heat flux data modified by the 0.76 factor given in Equation 58 (§ 3.1.2). Hereafter, in 
this study, these data were taken to be the critical heat flux for ignition for each specie in 


the two grain orientations. 
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Table 11. Final critical heat fluxes for ignition used in this study obtained by the analysis 
of the time to ignition data. 


For Redwood, Janssens [33] quotes a critical heat flux of 14 kW/m’ whilst Tran & White 
[34] give a value of 12.42 kW/m”. In this study, the critical heat flux for ignition of the 
Redwood along the grain was found to be 11.74 kW/m”. Janssens [33] gives a critical 
heat for ignition for Douglas fir of 13 kW/m? compared to 12.19 kW/m’ in this study for 
the along the grain samples. Tran & White [34] quote a critical heat flux for ignition of 
10.53 kW/m? for Red oak compared to the value of 8.20 kW/m’ for the along the grain 


orientation samples tested in this study. 
An overall comparison of the critical heat flux values obtained in this study compared 


with literature values show slightly lower values for the along grain orientation and 


significantly lower values for the across grain orientations. 
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Figure 30. Comparison of calculated and measured critical heat fluxes with data given by 
Janssens and Tran & White. 


The differences in the values may be partly explained by natural variation in the wood 
species but also by the fact that Tran & White's tests were conducted in the OSU and that 
Janssens [33] tested his samples in the Cone Calorimeter in the vertical orientation. 
However, in the study by Atreya et al. [21] it was found that the critical heat flux only 
varied by about 10% between horizontal and vertical samples and the critical heat flux 
was greater in the vertical case. Thus we might expect Janssens critical heat flux data to 


be somewhat less if his samples had been tested horizontally. 


Finally, the critical heat fluxes obtained in the experiments (Table 9) were compared with 
those derived from the analysis of the ignition data (Table 11). Figure 31 shows a plot of 
the experimental critical heat fluxes and derived critical heat fluxes. It is clear that 
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although the two do not match there appears to be a relationship between the two as 


shown by the linear regression fit. The difference between the derived and experimental 


critical heat fluxes may be as a result of neglecting the glowing ignition data for the low 


heat fluxes. However, since in only two cases were the no ignition incident heat flux 


found (Table 9), this comparison is not conclusive and further tests should be conducted 


to verify these findings. 
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Figure 31. Comparison of measured and calculated critical heat fluxes. 


3.2.3 Ignition temperature 


a) Average ignition temperature 


By obtaining the critical heat flux for ignition for each species in the along and across 


orientations, Equation 35 can be used to solve for the average ignition temperature. Since 


Equation 35 assumes a long ignition time, i.e. a low heat flux, the value for the 
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convective heat transfer to the surface of the sample in the Cone Calorimeter was taken 


as that at the critical heat flux using a correlation provided by Spearpoint [45] such that 


h, =—2.04g/+19.58 
Equation 67 
Equation 35 was then solved numerically by an iterative process for the critical heat flux 


to obtain a value for the average ignition temperature. Table 12 shows the calculated 


ignition temperatures obtained for the various species in the across and along grain 


configurations. 


Table 12. Calculated average ignition temperatures. 


The average ignition temperatures for Redwood and Douglas fir (softwoods) are 
generally greater than those for Red oak and Maple (hardwoods) in the two grain 


orientations which agrees with the findings discussed in § 1.2.2. 


83 


The average ignition temperatures obtained in this study were compared with data quoted 
in the literature. Tran & White [34] measured the ignition temperature of their samples 
with a thermocouple on the exposed surface of the samples. They quote an average 
ignition temperature for Redwood as 364 °C. Janssens [33] gives an average ignition 
temperature from Redwood as 363 °C. Dietenberger [36] gives ignition temperatures of 
353 °C in the Cone Calorimeter and values between 290 °C and 356 °C (depending on 
the moisture content of the samples) in the LIFT [11]. All of these values compare 
reasonably well with the average temperature calculated in this study for the along grain 
oriented Redwood with the value given here being slightly above those quoted by the 


other researchers. 


Janssens [33] quotes an ignition temperature of 350 °C for Douglas fir which is lower 


than the temperature of 384 °C calculated in this study for the along grain orientation. 


Tran & White [34] obtained an ignition temperature of 315 °C for Red oak and Atreya er 
al. [18] quotes 365 °C. Both of these values are greater than the ignition temperatures 
obtained in this study for the along and across grain orientations. The data from the 
literature and this study demonstrate that there is a fair degree of variability in the 
ignition temperatures of wood. As discussed earlier, there are many factors that influence 
the ignition properties of wood. However, the average ignition temperatures obtained in 


this study are comparable with the data quoted by other researchers and an average 
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ignition temperature of somewhere between 300 °C and 380 °C for along grain oriented 


wood Is typical. 


b) Ignition temperature as a function of incident heat flux 


By rearranging Equation 33 we obtain 


t, (gi) 


T, =T,+ |= 
i C,,kpc 


Equation 68 


Using the measured times to ignition and the apparent thermal inertia obtained in 3.2.4, 
the ignition temperature at a given incident heat flux can be calculated with Equation 68. 
Since Cig also includes Tig, Equation 68 cannot be solved analytically but has to be solved 


iteratively. 


In the study by Hopkins [30], thermocouples were also located on the exposed surface of 
the samples tested so as to obtain the ignition temperatures at given incident heat fluxes. 
Figure 32 compares these calculated ignition temperatures obtained from Equation 68 for 
Redwood with those measured by Hopkins [30] and quoted by Tran & White [34] and 


Janssens [33]. 
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Figure 32. Predicted ignition temperatures for given incident heat fluxes for Redwood 


The calculated ignition temperatures, particularly for the along grain configuration, 
compare well with the other literature data at heat fluxes above around 20 kW/m’. The 
ignition temperature quoted by Hopkins at 21 kW/m? is greater than those found 
elsewhere. Below 20 kW/m’ the calculated ignition temperatures show a downward trend 


with a limiting value of around 200 °C for the across grain configuration. 


Similarly, for Douglas fir, the calculated ignition temperatures at given incident heat 
fluxes are shown in Figure 33. As for Redwood, the édiciilated temperatures compare 
well with the data quoted by Janssens [33] at heat fluxes above 20 kW/m”. Again, below 
20 kW/m’ the temperatures decrease to values around 200 °C for the across grain 


configuration. 
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Figure 33. Predicted ignition temperatures for given incident heat fluxes for Douglas fir. 
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Figure 34. Predicted ignition temperatures for given incident heat fluxes for Red oak. 
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Figure 35. Predicted ignition temperatures for given incident heat fluxes for Maple. 


The calculated ignition temperatures for Red oak and Maple at given incident heat fluxes 


are shown in Figure 34 and Figure 35. 


In all four cases the ignition temperatures are almost constant at incident heat fluxes 
above around 20 kW/m’. Simms [15] quotes work by Bamford, Crank & Malan in which 
it is suggested that at high incident fluxes the energy required for surface ignition 
appeared to tend to a constant value. The results from this study agree with these 
findings. Except for Red oak, below 20 kW/m’, the ignition temperatures fall to values 
lower than the constant values found above 20 kW/m”. These results appear to confirm 
that the ignition mechanism is different at low incident heat fluxes as observed in the 
experiments (§ 1.4.3). A similar decrease in the measured ignition temperature of PMMA 


with decreasing incident heat flux was obtained by Rhodes & Quintiere [46] 
88 


The fact that the ignition temperature falls as the incident heat flux is reduced initially 
appears to conflict with Atreya et al. [21] in which they found that the ignition 
temperature rises as the incident heat flux decreases. However on close examination of 
their data (for Mahogany) it was found that the minimum incident heat flux used in their 
experiments was ~18 kW/m7”. The data obtained in this study for Douglas fir, Redwood 
and Maple shows that around this same flux region the ignition temperatures also shows a 
slight rise (Figure 32 and Figure 33) before decreasing again as the incident flux is 
further reduced. The ignition temperatures obtained by Hopkins [30] for Redwood also 


shows a rise at 21 kW/m? compared with 30 kW/m? and 42 kW/m’. 


It is interesting to note that the ignition temperatures of around 200 °C in this study at 
low incident heat fluxes is similar to those given for self-ignition temperatures of wood. 
For example, Cholin quotes [3] self-ignition temperatures of small samples of wood are 


quoted as being between 192 °C and 220 °C. 


3.2.4 Thermal inertia 


The apparent thermal inertia can also be obtained from the slope of the best-fit line of the 
plot of 1/ Ae against incident heat flux. From Equation 31, at ‘high’ heat fluxes (as 


defined in § 3.2.2), 
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Equation 69 
thus 


Table 13 shows the apparent thermal inertia values calculated for the various wood 


species tested using the average ignition temperatures from 3.2.3. 


Table 13. Calculated apparent thermal inertia for wood species tested. 


Equation 14 states that 
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Equation 70 


where the thermal conductivity k varies by a factor of around 2.1 due to grain orientation. 
§ 2.4 shows that k is also a function of temperature (and moisture), thus Equation 70 can 


be rewritten as 
1~([fk(T)P 


Equation 71 
Consider a typical along grain wood with f= 1 and normalised thermal inertia and 


thermal conductivity (Tig) of 1 each. For the equivalent across grain orientation f= 2.1 


and assuming the ignition temperature is twice that for the along grain case then, from 
: 0.13 
Figure 13, the normalised A(7jz) may be typically 6 ~ 1.01. Thus the normalised 


thermal inertia would be from Equation 71, (2.1 x 1.01)’ = 5.2 times that of the along 


grain case. 


The data obtained in this study shows increased values for the thermal inertia of the 
across grain orientations although the ratios between the across and along grains vary 
considerably (Table 14). Only Douglas fir exhibits a ratio approximately equivalent to the 


5.2 suggested above. 
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Table 14. Ratio of across to along grain thermal inertia. 


The along grain orientation values for the thermal inertia obtained in this study are 
greater than those quoted by Janssens [33]. For Douglas fir and Redwood, Janssens gives 
0.154 kJ’.m“*.K™.s" and 0.138 kJ’.m“.K”.s" respectively. Tran & White [34] obtained 
the apparent thermal inertia of Redwood and Red oak as 0.073 kJ?.m~.K”.s and 

0.360 kJ’.m~.K™.s" respectively. Again there is a noticeable difference between these 
data and those obtained in this study. Tran & White’s data for Redwood is also only 
approximately one half that obtained by Janssens. As for the critical heat flux and 
ignition temperature comparisons, these differences may be due to several factors; 
previous samples were tested oven dry whereas the samples tested in this study included 
higher levels of moisture; previous samples were tested in the vertical orientation in 


different apparatuses; the natural variation of wood. 


3.3 Dimensionless ignition analysis 
The time to ignition against incident heat flux data can be plotted in a dimensionless 
form. From Equation 69 we can define dimensionless irradiance and dimensionless time 


as 
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lay" Gok 
= and 7, = og respectively. 
B, Ver . (Tr, -T,)° kpe 


Rearranging Equation 33 we obtain 


if =—C, keel, -7,) 


1S 


Equation 72 
dividing through by q”” 


m2 
qi 1 2 
—= C.kpc\T., —T, 
ae he lig sd ( ai ) 


Equation 73 
thus, substituting the dimensionless forms 


Or 


Equation 74 
A dimensionless plot of all the ignition data is shown in Figure 36. The plot also shows 


the theoretical curves with Cj, having either the 4/3 or 27/4 factors. The plot shows that 


scatter of the data is within the bounds of either the 4/3 or 27/4 factor used in the theory. 
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Figure 36. Dimensionless ignition plot on linear scales for all species tested showing 
comparison between measured ignition times and theoretical values. 


The data shown in Figure 36 was plotted on log scales to show the low incident heat flux 


data more clearly (Figure 37). It can be seen that the experimental data at low heat fluxes 


does not match the theory. The data does not curve as sharply to Vp =1 as the integral 
ig 

model solution suggests. This discrepancy between the data and theory is a result of the 

localised ignition mechanism observed in the experiments as described in § 1.4.3. In the 


theory we only account for the external heat flux and not any additional energy from the 


glowing process. 
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Figure 37. Dimensionless ignition plot on logarithmic scales for all species tested 
showing comparison between measured ignition times and theoretical 
values. 


From Equation 74 the gradient of the dimensionless plot gives 
fd 
eof 1) 
gradient 


By plotting a best-fit line through the data shown in Figure 36, a gradient of 1.21 is 


Equation 75 


obtained. Thus, from Equation 75, Ci, is found to be 0.68. The value for Cj, compares 
well to the 0.62 quoted by Abu-zaid & Atreya [20] and the gradient of 1.21 is close to the 


4/3 (= 1.333) value predicted by the integral solution. 
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3.4 Transitional critical heat flux 

As noted in § 3.2.1, the difference between times to ignition for the along and across 
grain samples for a given wood exhibited a distinct pattern. The across grain orientation 
samples ignited sooner at the lower incident heat fluxes (~20 kW/m’ or below). Yet at the 


high heat fluxes the along grain samples ignited quicker. This pattern can also be seen in 
the plots for 1/ aie against incident heat flux where the linear regression and theoretical 


lines intersect at some ‘transitional’ incident heat flux (i.e. the heat flux at which the 


along and across grain ignition times are the same). 


To find this transitional heat flux from the plots of 1/ lie against incident heat flux we 


note that, 


1/ Nips 1, = gradient, .g; , + intercept, 


Equation 76 
1/ tie. = gradient, gi , + intercept, 
Equation 77 
Thus, at the transitional heat flux, 1/ aes L= / Af x and qd, =4;, =4; and so 


substituting and rearranging Equation 76 and Equation 77 


., _ Intercept, — intercept, 
me A 


gradient, — gradient y 
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The transitional critical heat fluxes were determined as 21.7 kW/m’ for Douglas fir, 

25.7 kW/m’ for Redwood, 24.0 kW/m” for Maple but only 13.75 kW/m’ for Red oak. It 
interesting to note that for Doreredr Maple and Redwood the values for the transitional 
heat flux are similar. These values also correspond to the suggested boundary between 


the ‘high’ and ‘low’ incident heat flux values given in § 3.2.2. 
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4. BURNING RATE 

4.1 Integral model theory 

4.1.1 Assumptions 

The ignition process described in the previous chapter determines the initial conditions 
for the burning rate model. The assumptions used in the ignition model are carried 
forward to the burning rate model with the additional assumptions specified. Thus, for the 


decomposition model (Figure 38) it is assumed that 


(a) The fuel decomposes to gaseous fuel (volatiles) and char in an infinitesimal 
pyrolysis front at a fixed vaporisation temperature, 

(b) The solid is infinitely thick, 

(c) The virgin wood is inert up to ignition and decomposition, 

(d) The char material is also inert, 

(e) The flame heat flux remains constant, 

(f) The density of the volatiles is much less than the density of the virgin wood 
and the char. 

(g) Material properties are constant over the range of temperatures considered, 

(h) The volatiles do not accumulate within the char layer but are produced and 


exit immediately. 


For simplicity in assumption (a), we assume that the vaporisation temperature and the 


ignition temperature are the same though in reality these differ as compared to a liquid 
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fuel's flash point and boiling point. Suuberg et a/. [22] indicate that chemical kinetics are 
important, but also two separate kinetic regimes that must be included. They also find 
that the assumption of an infinitesimal pyrolysis zone is valid for an incident heat flux 


greater than 40 kW/m”. 


Assumption (c) tells us that it is presumed that there 1s no mass loss prior to ignition as a 
result of the pyrolysis of the wood due to the incident heat flux. The mass loss only 


occurs after ignition due to the flame heat flux. 


Assumption (d) implies that the exposed char layer is not affected by the incident heat | 
flux. However, as discussed in § 1.2.2, the char layer will crack due to pressure gradients 
within the material. The char will also undergo oxidation that will contribute to the mass 


loss rate of the sample and will reduce the dimensions of the char layer. 


a 


pyrolysis 
front 


| virgin wood 


Figure 38. Schematic of model for decomposition. 


4.1.2 Flame boundary condition 

The heat flux from the flame by radiation and convection can be intimately coupled to the 
burning rate of the mass supply rate of the fuel. In the Cone Calorimeter it is assumed 
here that the flame radiant flux to its base is constant for a given anteial provided the 


flame is tall i.e. its height H is greater than its diameter D (Figure 39). 


100 


Figure 39. Assumed cylindrical flame. 


The emissivity of a cylindrical flame ¢ can be given as 


—~1_, “sn 


Equation 78 


where x is the absorption coefficient which depends on the fuel and Lz is the mean beam 


length [46]. Thus, as Ly increases, & approaches a constant value. 
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Figure 40. Relationship between flame height, diameter and mean beam length (from 
reference 47). 


Meee L : 
From Seigel & Howell [47] and shown in Figure 40, an — 0.81 for <= >3. Since, as 
discussed by Rhodes & Quintiere [46], 9% rq =E_OTq ANd Gh com ¥ Constant , it follows 
that for = > 3, La becomes constant, thus from Equation 78, & becomes constant and 


hence g/,,.4 is constant. As an approximation we assume that the flame heat flux q% ,.q 1S 


always constant even during the growth or decay of the flame though it 1s recognised that 
the flame heat flux is coupled to the fuel supply rate which affects the flame height, 


thickness and surface convection. 


Just after ignition the net surface heat flux is the sum of the external incident flux and the 


flame heat flux minus radiation losses 
qt(7,)=4r+9,-o(T' -T), x=0 


Equation 79 
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The absorption of the incident heat flux due to the presence of the flame is not 
considered. Rhodes and Quintiere [46] showed that the flame from a sample of black 


PMMA tested in the Cone Calorimeter is almost transparent (>90%) to the incident flux. 


Thus just before ignition when t= t,,, 


from Equation 16 


q' =4" =q[-o(t{ -T;)-n,(, - 1) 


Equation 80 
and when ¢ =7;, , from Equation 79 
a = 4(0,)=44+4)-o(t - 73) 
Equation 81 


Hence there is a step-discontinuous change in the surface heat flux at ignition. 


4.1.3 Conservation equations 
The decomposition of the charring material is split into three regions; the char layer, the 
pyrolysis front and the virgin material (Figure 41). Each zone is considered in control 


volume form. 
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Figure 41. Burning rate model control volumes. 


a) Char layer 
The char matrix contains solid char and fuel gases. These are described in terms of bulk 


properties based on the volume V of the char-gas matrix such that 


m, =V.p, 
Equation 82 
where m, and pz are the mass and density of the gas respectively, and 
m, =V.p, 
Equation 83 
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where m¢ and /¢ are the mass and density of the char respectively, and 
m=m, +m, 
thus 


P=P,t Py 


Using conservation of mass on the char layer control volume we obtain 


“ As pas— Ay -i)+ A p,[ vp -1)- Ze} 


dé 
+A, ol vl) “eal i)=0 


‘Substituting Equation 84 into Equation 85 and dividing through by Ao 


#e do, do, 
Dye he) aS Ps TNS Dd ERT = 


thus, cancelling through we obtain 


el ie, 


Consider conservation of energy for the char layer 
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Equation 84 


Equation 85 


Equation 86 


Equation 87 


aes u,+ Py u, de — —m" h ml )(-i) 
ie (r.) - tif Selb ) 
+eph{t,) 0- Ze) (+i) 
+p, h,(T,)[0- o}(-#) 


a t. dé, .. 
=q"(T,)-45 ~ PUT, tis 


ad 
dt 


Equation 88 
where qj is the heat flux from the char to the pyrolysis front and from Equation 79, 


q. (7.) is the net surface heat flux at the temperature 7,. Rearranging we obtain 


da 
dé 
+p, hg(C,)— py hy(T,)+ 5m 
=G3(T,)- 4% 
Equation 89 
Given that 
h=u+l 
Equation 90 
we can express 
T, 
(1) =u, + PE) 
Ps 
Equation 91 


and 
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n(T,) =u,(r,)+ 22) 
Ps 


Equation 92 
We can use Equation 91 and Equation 92 to write from Equation 89 


ae h(E) ea) o.) =|)» PED], )» PED] 


=-p,u,(T,)- Pp, sts p,t,{T,)= Py. PE). of) 
=-p,u,(T,)— p,u,(Z,)— p(Z,) 


Thus substituting back into Equation 89 


d # m 
dt is (0,1, + pyu,)dx + ri" h,(T,)— ri” (T,) 


th p,u,()- eyup(tIS* = a2(t,)-45 - pl) 


Equation 93 
From thermodynamics, let u, =c,,T, u,=c,T and take h, =c, , T . Assuming that 


Cy =C,, =C, then Equation 93 can be written as 


a 


al * (p, ¢,T + Pj C,T)dx + mn"c, T, —m"c, T, 


OSs ae. 
+- Recel, = 447, * = 4 — 4% 
Equation 94 
Since 


d 


a Cgc, 7, — P57, Jax 


dé. 
- P. eel, — 2467, | = 


we can rearrange and write Equation 94 as 
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v 


at (lo, ¢, T+ pe, T— (ep ed ~ 1,6,T,)dx+ m' c, (T. -T,) 


=q"(T )-—q" + D(T aor 
g’( 5) qi, + Pl y) dat 


Equation 95 
d 5, 
dt Mle. e+ P4e4} ~(p, ¢, + P4e,)T, \ax +m", (T, -T,) 
- - dé 
a qi (T,)- q+ p(T, )—* 
Equation 96 
5 
dr - yt ta dé 
(0, Cet Py “)s |v -7,)ax + mC, (7, -T,)= g'(T,)- 45 + p(T, )—* 
0 
Equation 97 


Typically Pz, << Pgs since Pz ~! kg/m’ and Cz *! kI/kg.K while Ps ~ 200kp/m° 
and Cy, * 3kJ/kg.K. In addition, 2, (7-T,)>> P(T,) since (T-T,)*100K and 
p~10° N/m’, thus 24% (T-T,)* 6x10" ky/m? or 6 x 10’ N/m”. Thus Equation 97 can be 


written as 


Equation 98 
b) Pyrolysis front 


Using conservation of mass on the pyrolysis front control volume we obtain 
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$I4 makers h, ie )-Seale +4, a -Seale / 


ce ES, 
=0 


if p.| -Se}t )=0 


Equation 99 
Thus, dividing through by Ao 
i dé, dé, do, ay 
eA aie (re ea 
do, 
Pw — Pg — Pe at = Pz 
Equation 100 
Substituting Equation 87 from the char layer analysis into Equation 100 
dé, 5 
(0, —P¢ “Ae ae 
Equation 101 
Equation 102 


Consider conservation of energy for the pyrolysis front 
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sla u, dx+ p,h, ) (-i)- eC i) 


do 


+0, hee ) 
+pahs(f} 0-2e [ei 
= (5-45 “wt SO} -)- ve | 
| Equation 103 


dé, 
a%s 4 h, (T,)+ p, h,(T,)—* 


dt t+ pyh, (T)*— ph TJ = 45 - ah 


Equation 104 


is dé sn 
ra’ h(T,)+—* [p, hy (T,) + Pg ly L,)~ Prob DN) = 45 — 4 


Equation 105 
Substituting Equation 101 into Equation 105 


do. dé re 


do “noon 
Bi pi) h,(T,)— 2g hg(T,)+ Pg he (T,) + Pp Me) bt oa de 


thus 
pel f(t) helt) (0) 0) |= asa 


Equation 106 
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Based on the mass of the virgin wood, we define the heat of vaporisation or heat of 


pyrolysis as 


a, = f(7,)-h.(,)}- 24 (0,)- (0) 


w 


Equation 107 
then Equation 106 becomes 
p Wee a A+ lan, = ge qi, 
Equation 108 


c) Virgin wood 


Consider conservation of energy for the virgin wood control volume 


& fo. u, dx +p, h,(T, yo-" os OG i)+ py hy (To Jo 2 At, 32) 49.) fi } () 


dé 0,+06 
mt, ~45-plr,) ot as rt) O22 


Equation 109 
Let 4“, =¢, 7 and assume negligible heat loss to the back i.e. q, =9 
54+, 

pice dé Ave +6,) 

ot pre, Th)d+ ph (Dp, kG) 

dt J dt 

Me dé, d\6,+6 
= 4+ lt) pf, ete) 

Equation 110 


Expanding for hy(7y) using Equation 90, then 
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27 pestle py [ut)+ AED] ep, ;)» ED] O28) 


= 41+ p(n) PT, ease 


Equation 111 


Since 4, =, T and the P(T) terms cancel then Equation 111 can be written as 


54+6, 
dé, d\6,+6 
— | p,c,T(x)de+ p, c, T,.—* yr rae oT ete a 
Wel 
dé, do, dé, 
T(x)dx oe 7 
a Pu ly T(x) + Py Cy Ty" — Py by Ty" = Pu Ov ToS" = Gh 
aay Piers 
; J eeecTelde+ a,c a <r, -T,)= Pye, — q oa q. 
Equation 112 
Let ¥ = ~x'+9, then Equation 112 becomes 
p,, 4 re.)-Tle+ PC 4 (7 -1, )=q" 
w ” dt ; 0 ww dt v 0 w 
Equation 113 


4.1.4 Temperature profiles 
Taken from the analysis by Anderson [29], temperature profiles for the char layer and 


virgin wood are selected that satisfy the specified boundary conditions. 
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a) Char layer 


The temperature profile in the char layer is assumed to be linear function with depth with 


boundary conditions such that 
(1) when x =0, T=T, 
(ii) when x=6,, T =T, 


Consider a profile such that 


x 
r-T,=(, -rj1-2| 


¢ 


When x=0 
0 
r-1.=¢,-T)1-2| 
Oy 
3 
thus 
i hey 


and boundary condition (i) is satisfied. When x = 6, 


(3) 
T-T,=(1,-7,}1-+ 


thus 
and boundary condition (11) is satisfied. 


113 


Equation 114 


b) Virgin wood 


The temperature profile in the virgin wood is assumed to be some function of depth with 


boundary conditions such that 

(i) when x'=0, T=T, 

(ii) when x'=6,, T=T, 

(iii) when x'=6,, ee =0 
Ox 


Consider a profile such that 


; 2 
x 
ToT (Te Ti 


Ww 


Equation 115 


thus 


hence boundary condition (i) is satisfied. When x'= 6,, 


2 
r-m=(,-7,|1-5) 


=0 
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thus 
T=f, 


hence boundary condition (ii) is satisfied. By expanding Equation 115 we obtain 
Qxete 
Liat RAE BE Merete 


Equation 116 
and differentiating with respect to x' 


hoy on ch aiettal aoe 


w Ww w Ww Ww 


oT 0%, _ AT, _, AT, (1), a7, (2x')_a% , ,O% (1 )_ ah ( 2x' 
& a& a a'l6,) alse?) ao a'\6,) a&\e& 


since x'=6,,, by substitution 


Oumecmectn of | 1 | eT, (20, \ oT, 07, (1 \_ at, ( 28, 
BO Gn ee OSG | Gall Bey Gal eo ir ex | 


' 12 ’ ! 12 
2¢-n)-2 [72-42% r, 2 7 Te) 


thus 


Gl sO] Cl eor, 


Bo ES eS 


ar _ of, 
Ova Ox. 


and since 7, is a constant then 


which satisfies boundary condition (111). 
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4.1.5 Governing equations 


Now that we have considered the conservation of mass, conservation of energy and 


temperature profiles for the char layer, pyrolysis front and virgin wood we need to 


develop our governing equations. 


Given that 


Differentiating Equation 114 


Clam iS ek, 
Fee n| 3) 


and substituting into Equation 117 


Differentiating Equation 116 


2 (,-7-2+2| 


116 


Equation 117 


Equation 118 


Equation 119 


Equation 120 


thus when *'= 0 


Equation 121 


Equation 122 
From the analysis of the char layer, Equation 98 states 


7) 
Peg \(T-T, dx + mc, (T, -T,)=95(T,)- 4% 
0 
and from the temperature profile in the char layer given in Equation 114, at time ¢ 


x 
Pde -{1-2] 


¢ 


thus 


fe- T,) dk = ja- nl -2}s 


Equation 123 
substituting Equation 123 and Equation 119 into Equation 98 we obtain 
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Poo d k ue T,) 


5 ath -T,)5, + m'c,(T, -T,) = 4" (T,)-+4+— 


04 


Equation 124 
From the analysis of the pyrolysis front Equation 102 states 
dé. 
=p Lae 
Be. mls 
substituting for the char fraction, % = ~ 
dé. 
m’ = p.(1—-¢)—* 
pull-4)— 
i _, A, 
(-¢) at 
Equation 125 
Equation 108 states 
dé 
—*AH, =q"-@" 
ae i = 5 Vy 
substituting Equation 108 with Equation 119 and Equation 122 
a4 an = K (7, SHE) 2k, (T, -T,) 
shay? | ph a 6, 6, 
Equation 126 
substituting Equation 125 gives 
tay ba) ke ad) 
(1 vs ¢) : om , 
Equation 127 


From the analysis of the virgin wood, Equation 113 states 
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pic. firle.s)-Tlas'+ p, Peet, -1)= 4 
dt : w at w v w 
from the temperature profile in the virgin wood given in Equation 115, at time ¢ and x’ 


; 2 
TET AT eT ee 
s=@.-1) 1-2] 


Ov. Ou »\2 
[I7@.)-T]&'= (7, 71-2) dr’ 
0 0 O, 


Equation 128 
substituting Equation 120, Equation 125 and Equation 128 into Equation 113 


(T.-T.)d5, mm" 


eee SE 
Pwlw 3 dt *T ay! v 9) 


2k,(Z, -T) 
re 


w 


dividing through by Pw“w (7, es 7) 


m 


lecraaet (<6) ann 


Equation 129 
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Adding Equation 124 and Equation 127 we obtain 


Poo a 


* = oy the —T,)5, + 0h"c, (T, -T,)+ => 


[aH] 


i ) 


ey BEAT) BD) oie 
=i) PMI Pe 
¢ ¢ w 

Poo d 


2k, (7, -T) 
2 tat O. 


GiB, +e, (7, 7) +7 [aH ]=45(0,)- 
Equation 130 
n” en(r\_ 2Ky(T, —Ty | 
ry tteg tT) + la =O 


w 


Pees 27,4 
Dadi di 


Equation 131 


If the time just after ignition is considered then t =t,," with T, = T, = Tig, m" = mj;, and 


Ow = Og. Thus, from Equation 131 


Psa(d. da i 7) 2kulf - 1) 
2 (Sr, Prt rites) +E plaed= G(T.) és 
=0 
hence 
m" _(i-¢) , 2k, (T,-T) 
ca, SH ale Oig 


Equation 132 
Assuming that at ignition 7, = Ti, 6,, =6,, and q', =q" then from Equation 122 


_ 2k (T=T) 
q” 


Equation 133 
Substituting for dj, in Equation 132 
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from the ignition analysis, Equation 16 can be used to approximate 4" such that 


q’ =a", )=gf-o(T! -T!)- AE, -T) 
and from Equation 81 
gf’ =a", )=0,+4'-off! 


Thus substituting into Equation 134 we obtain 


ni, = 9M) los + @t-oltt —n)-Gr-o$ 1), -7)} 
hence 
it (Dion ace, -T)| 


AH, 


It is interesting to note that Equation 135 suggests that the mass loss at ignition does not 


explicitly depend on the incident heat flux but only on the flame heat flux Gn . This fits 


with the assumption (c) given in § 4.1.1. However, it will be shown later that the char 


fraction ? and the convective heat transfer coffeicient h- can be considered to be 


functions of the incident heat flux. 
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4.1.6 Dimensionless analysis 
Finally we make the governing equations dimensionless by defining the following 


dimensionless variables specified by Anderson [29] 


m” L m;, L 


Vie ee ——_2__ 
(1—¢)q°(Z,) “¢ ei 
KSeae yee Nap tas 
yen th: ieee a ab 
T, v y 
t 
ra; Aja 1a) 
s Pw 
2k. L 
Oe ee 
c, @(T,) 


plus the dimensional variable 
L=AH, +c,(T,-T) 


where L is the heat of gasification per unit mass of virgin wood. 


From Equation 130 


Pale d 2k{T,—Ty) 
7 rer -T,)5,+m"c, (7, =e ire aH l= q'(T,)- Be 


Substituting for the dimensionless variables M, 7, 6;, 09, A and Ag 


pete M (1-9) 4» 3 
2 rEvaeneeL 7 aT Jeg 7.8, -1) 


CD arr) ME = g(r.) EEA) 


1972 


re) 
Multiply through by LT to obtain 


8, Ps% ky 5, M(-¢) 
6-1 a 
k, (ae 2 Le Cy }, at sale 7, (6, ee har iy Mi(rJe, T, (0, 1) 


+O MAM jo(p) 8. gr(r) 5: 2k, T,(1-) 
ae L kT ree Xo, 
heh 
: Ye) = 
Since 9; c, 9° (T,) 
1( Pc | d 2k,L M(i-¢) 
@.—-1)A Tr and 
3 Get gle ade ty ae Fil 
Qin le wie MAH pA od EO! Tod lions 2(1-4,) 
7a Cae *9"(T,)=——.— #"(7,)- 2 
c, 9" (T,) kT, i qn i) : "(T.) kT. qr ( i) A 


3 2} “fe -Na,]e2{ Se, 0a) 


ae 
fe 2 g(t.) 20-6) 


ut CAVE) VMS 


wv 


Thus substituting for ? = p — and rearranging, 


w 


Equation 136. 
From Equation 127 


wt gry HoT) _ 246-7) 
(1-¢) : O. é, 


substituting for M, 44, A 8 and % 
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M (\-4)4'(7,) AH 
jk (i- (-¢)~ Lae 


and then substituting for % 


MG@i(T,)AH, _, co G(Z) 1 rm 1) are 6 4.(Z,) | welll 6,)T, 
2k, L 


ib Dir yas ae 


Rp maTNTycy alk ao to) 


g'(T.) AH, pb Ny 


From Equation 129 


k, 
doom tT sugary 


1 
Bat i pl=s)) 6, 


substituting for M, 4 and 7 


1 a(ad,) , MQ-9)a(Z) ee ea 
3 { 26: L “Pull-9) Pet, AS, 
a 
ldda M@(T,)_ 2k 1 
Sar oul DAL tel “AO, 
1 dd , 6, Puy MGT) _ 2ky Py ly 1 
3 dt ke peL Dic eke a 


A 


1dA 2k,L c,M@'(T,)_ 2 


we G(T.) Seeker 


thus 


A 
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is —— (0, -1)T, -2k Fan 6,)T, 


| 


Equation 137 


Equation 138 
a) At ignition, from Equation 134 
.n  (l-@)ypew ca 
ma, = ral dar(r.)-4"] 


v 


Substituting for the dimensionless variable Mig 


EU Pn) 4") 


®” (-¢)q" AH, 


Equation 139 
From Equation 27 


Paty ia B 
Oig on _n 
Since Ai, = So” at ignition ¢ = fig and 9 (7,)= q(T), thus 
2-8 
3 cl eerie t, 
Rac) Ane 


ig ~ 2k, L 


c, 7" (T,) 


Equation 140 
From Equation 28 


4 ope (7, -7,) 
i, =k, Pw [w : 
~ 3 ; 2 Ps 
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substituting into Equation 140 and cancelling common terms 


ee 


qr(7,) 


lea 


q ig 


Cw (7, -T,) 
L 


Equation 14] 


b) Consider times just after ignition (i.e. small-time). From Equation 125 


mm" ua 


(-9) ” 


Substituting for the dimensionless variables M, 7 and 4¢ 
MGT.) _ _(A,6,) 
ip df) 6 if 
M@i(T,)__ky As 


L  c,o, at 


Substituting for % 


Lik 2 Sl) 


yu we 
vear 


At time “ig, char depth 4¢= 0 and at time 7, char depth is A¢, Integrating Equation 142 


Equation 142 
is A, 1 
[Mdr= [—aa, 
Tig 0 2 
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M(c-r,,)= ay 


since M@ ~ M,, at Tz then 


A, ~2M,, (c—r, ) 


From Equation 136, with no char 4, = 9 and assuming T, =T, thus 9, =1and 


Equation ]43 
q3(T,)> 4.(7,) 
we oy 4l6 -1)A + (1-¢) Lee (@ ~1)+ AH, M= L q'(7,)_(i-4) 
denteoablde {- hy = coe) i c, T, cul) g: (1s) eA 
Ne pe! ee —S>— 
=0 =0 4 | 
peimescsT=0,) 
Alle AH Ba A 


Equation 144 
E 
: 6, 2-2 
Since % e then 
Ty 
¢, T, (1- @,)=«,7,| 1-3 |=, (%, -) 
substituting back into Equation 144 
ryt ine AAU cade ah 
AH, AH, A 
Equation 145 
alternatively 
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Equation 146 


substituting for M using Equation 145 and rearranging 


ae pcm 4-4-0 


3dr | AH, AH, A A 
Equation 147 
Since, from our definition of L 
c, (7, -T))=L—AH, 
Equation 147 becomes 
1da |_L _(L-AH,)1_1)_, 
3 dt AH, AH Aga A 
1 dA L Ee le Al Pelee 
——+2 - —+—.—-— |=0 
3 dt AH, AH, A AH, A A 
=0 
dA __ _L -2| 
at AH, A 
ds _{_L 4] 
dat AH, A 
Equation 148 


Let ¥Y=1-A 
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a a 


v 


y AH, 
Equation 149 
by integration of Equation 149 
jig Z 
ydt dt AH, 
VY 
In| — |- —6 TT, 
[¥) -o- mg) Eo r—re) 
on = on - A)- 6 Je eke ) 
ig v 
Equation 150 


Since at early times 4— Aig is small, 7 — Tig is small and e* =1+X for small x, Equation 


150 can be approximated as 


by rearranging 


L 
ba1-(-a,)+(-ae)0-a,)6f 2 


=0 


Equation 151 
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c) Finally consider times much later than the ignition time (i.e. long-time) 


dA 
Let the char depth rate approach zero, i — 0 the surface temperature approach a 


dé. ate 
constant, "Ey —>0 , and the burning rate tend to zero, M > 0. From Equation 138 


3 dt 
Bidtack 


By integration between ignition and time 7 


[aaa — [eae 
A, 


ig Tig 


~ 2 — 
Aw pi +12(r Tip 


Equation 152 
From Equation 136 


14%] 2(e.-a,)+|0-0) Se} - 2 Ju OE} 


From Equation 79 
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G(T.) = af +d -o(T! -T)) 
thus by substituting, rearranging and dividing by ig 


a(t) -T - Et Oak a) 
7 Te (per A 


v v 


As A> ©then 


OWT -Ty) Gn +4 
as ae 


v 


substituting for % and % 


Equating M’s in Equation 137 and Equation 142 


1dA, _@t(7,) Tc, (22-9) 


QAP gt AHA hen meas A 


assuming /; =, thus g(7,)> 4(7,) 


sabe Ree (Fe ne) 8) 


2 dt AH,\\k,J2 A, A 
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Equation 153 


further assume % is constant, the thermal penetration depth approaches infinity, A>, 


as time tends to infinity, 7 > © and finally at time Tig the char depth 4¢= 0 


Equation 154 
From Equation 142 
_1 a, 
2 dt 
substituting Equation 154 
ina Cals 
M =~ ——| ,[2| *— |6(@, -1)\r-7, 
el Ape 
thus 
Equation 155 
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Thus the dimensionless groups of relevance are as follows - 


Ignition: 
M.= em ya q- (Equation 139) 
IAT ang. (7) 
A,, = caller) chi (r,,) (Equation 141) 
L q; Lie 


Small-time burning rate solutions: 


M= (= [ - aaa ' (Equation 146) 
AzxA, 4% Jan en) (Equation 151) 
A, ~2M,, (c - Ti) (Equation 143) 


Long-time burning rate solutions: 


2c, Haas be T— Ti 


A, * Sun Bio AT Man, (from Equation 154) 
oc, \T. -T, : 

M~ ioe eee from Equation 155) 
SAH, (7 —Tig ( z 

Ax JA, +12|z —7,, . (Equation 152) 
7 aM 1/4 

T, = [a4] (from Equation 153) 

o 
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As discussed in § 2.6, we can replace the ratio of the thermal conductivity of char and the 


thermal conductivity of wood 7A with the char fraction @. 


4.1.7 Solutions of the integral model 

The general form of the dimensionless equations is shown in Figure 42 and the equivalent 
dimensional variables in Figure 43. The small-time burning rate solution is equivalent to 
what would be obtained by a non-charring material as detailed by Hopkins [30]. At 
ignition the small-time burning rate starts at a given point after which it increases to a 
maximum steady-state value, where it remains. A real burning material would of course 


eventually consume all its fuel and the burning rate would drop to zero. 


The long-time solution effectively starts from an infinite burning rate from where it 
decays towards zero. There is a transition point at which the small-time and long-time 
solutions cross (marked by the diamond-shaped symbol in Figure 42). This point marks 
the peak burning rate and it is less than the steady-state non-charring material burning 
rate. Thus for a charring material the burning rate follows the small-time solution up to 


the transition point (the peak burning rate) and thereafter follows the long-time solution. 
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1.6 


small-time RR ie Sage Mossy 
| solution eRe cae 
1.4 ig RS 3 
| eee 
ea: 
= }e2 | Ef 
= l Hes 
a LO om 
ie \ pai 
oe y 
= AQ Mpa 
Ss 0.8 / \ 
2 / 
Es M, \ 
S 0.6 ~ 
5 pix 
£ ~ 
B 04 a 
Get et long-time 


0.00 0.10 0.20 0.30 0.40 0.50 0.60 0.70 0.80 0.90 1.00 
Dimeasionless time, t - t,, [-] 


Figure 42. Dimensionless burning rate solutions. 


m 
steady 
0.050 
1 | qaeens ee oP a ee 
Saito ee eh Nr eee 

0.045 | | Pies a= 

| Ne 
0.040 | ae 
0.035 | | se 

i Me 


0.030 ae 


0.025 9 ™" peak 


Burning rate, m"' [kg/s.m?| 


0.020 m’. 
\ 
0.015 \ 
| \ 
0.010 > Sion rai 
Fi Se solution 
0.005 
0.000 
0 t. 60 120 180 240 300 360 
ft Time, t {s] 


Figure 43. Dimensional burning rate solutions. 
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Other parameters such as the char depth and thermal penetration depth exhibit a similar 
behaviour. For the thermal penetration depth we can also include the period prior to 
ignition. Figure 44 shows the thermal penetration in dimensionless form and Figure 45 
shows the thermal penetration in dimensional form. Figure 45 has been re-plotted in 
Figure 46 to show more clearly the transition between the pre-ignition, small-time and 
long-time solutions by removing the non-applicable portions of the curves. The 
dimensionless and dimensional char depth solutions are shown in Figure 47 and Figure 


48. 


Dimensionless thermal penetration depth, A [-] 


0.20 0.00 0.20 0.40 0.60 0.80 1.00 1.20 1.40 1.60 1.80 2.00 
Dimensionless time, t - t,, [-] 


Figure 44. Dimensionless thermal penetration depth solutions. 
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Thermal penetration depth, 5 [m] 


0.025 


small-time oie 
solution ee one 
NS, toni 
Pa is _— 
0.020 long-time eet 
solution LH 


0.015 


0.010 


0.005 


0 30 60 90 120 150 180 210 240 
Time, t [s] 


Figure 45. Dimensional thermal penetration depth solutions. 


Thermal penetration depth, 5 [m] 


0.025 . 
long-time ora Oe 


_ 


solution SS jah ee 
-_ 
= 
—_ 


0.020 
0.015 
0:010 


0.005 


0.000 
0 30 60 90 120 150 180 210 240 - 


Time, t [s] 


Figure 46. Dimensional thermal penetration depth solutions re-plotted for clarity. 


137 


Dimensionless char depth, A, [-] 


0.00 0.10 0.20 0.30 0.40 0.50 0.60 0.70 0.80 0.90 1.00 
Dimensionless time, t - Te [-] 


Figure 47. Dimensionless char depth solutions. 


0.008 , 
0.007 os 
| ,/’ small-time 
0.006 x solution 
/ 
= 0.005 Ws 
¢ é are 
a 0.004 + y Me ihc 
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a Be et solution 
5 Seog k Sie ao 
/ eS has 
0.002 | ae 
i —_ 
eg 
0.001 | by be 
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Figure 48. Dimensional char depth solutions. 
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We can now examine these dimensionless groups to investigate the form of the small and 
long time solutions. Each dimensionless group is a relatively complex function of many 
variables. In order to examine how each variable contributes to the behaviour of the 
model, we will select a set of typical values for these variables and then vary one at a 
time. The values were chosen on the basis of properties measured for a typical wood (§ 2) 
and those obtained by the analysis of the experimental data (§ 4.6.1) and are presented in 


Table 15. 
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50,000 
20,000 


15,500 
60,297 


36,949 
0.310 


0.00000021 
1,500,000 


803,207 


Sia 
iss 


Sis 


12,000,000 


| 15500,000 


5.67E-08 


20 
200,000 


Table 15. Typical values used for variable analysis. 
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The effect on the predicted burning rate and the predicted char depth are considered. In 
each case, the small-time and long-time solutions have been varied above and below the 
typical value. Plots using the dimensionless variables are presented, however, these can 
be difficult to relate to real world measurements and thus the plots are also shown in 
terms of dimensional values. Common x-axis and y-axis scales have been used for each 
plot to allow for direct comparisons to be made between the relative effect of each 
variable. In addition, each plot shows the transition point (marked by the diamond-shaped 


symbol) between the small-time and long-time solutions. 


These dimensionless plots present an idealised view of the effect of each variable. As 
already discussed, wood is a complex material in which its properties vary between 
species and between individual samples. Therefore it may not be so obvious which 
properties are having the greatest influence on the burning rate of a particular sample 


when tested experimentally. 


a) Incident heat flux, 7 

Plots for the variation in the incident heat flux (Figure 49 to Figure 52) are shown at the 
values used in the main experimental test series, namely 25 kW/m’, 35 kW/m’, 50 kW/m” 
and 75 kW/m”. Since the char fraction was found to be a function of the incident heat flux 


(§ 4.5.1), this varies as labelled. 
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The analysis shows that as the incident heat flux is increased so is the burning rate and 


the small-time solution is more sensitive to the influence of the incident flux. The 


analysis agrees with the findings of the work conducted at the Joint Fire Research 


Organisation (JFRO) [48] which also concluded that the burning rate increases with 


incident flux. 


The char depth solutions show that at low incident heat fluxes the small-time solution 


predominates over the long-time solution. In Figure 51 and Figure 52 the long-time 


solution for g” = 25 kW/m’ does not appear in the time scale shown. 
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Figure 49. Effect of incident heat flux on burning rate, dimensionless plot. 
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Figure 50. Effect of incident heat flux on burning rate, dimensional plot. 
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Figure 51. Effect of incident heat flux on char depth, dimensionless plot. 
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Figure 52. Effect of incident heat flux on char depth, dimensional plot. 


b) Flame heat flux, q7, 


The flame heat flux was taken to be values of 15 kW/m?, 20 kWim?, 25 kW/m? and 


30 kW/m’ and the results are shown in Figure 53 to Figure 56. 
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Figure 53. Effect of flame heat flux on burning rate, dimensionless plot. 
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Figure 54. Effect of flame heat flux on burning rate, dimensional plot. 
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Figure 55. Effect of flame heat flux on char depth, dimensionless plot. 
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Figure 56. Effect of flame heat flux on char depth, dimensional plot. 
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C) Critical heat flux, q’ 

Plots for the variation in the critical heat flux (Figure 57 to Figure 60) are shown at 
values of 15.5 kW/m, 14 kW/m’, 12.5 kW/m’ and 11 kW/m’. As discussed in § 3.2.3, 
the ignition temperature of the material is a function of the critical heat flux and thus the 
calculated values for the ignition temperature are indicated on each plot. As with the 


incident heat flux, the char fraction also varies with critical heat flux. 
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Figure 57. Effect of critical heat flux on burning rate, dimensionless plot. 
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Figure 58. Effect of critical heat flux on burning rate, dimensional plot. 
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Figure 59. Effect of critical heat flux on char depth, dimensionless plot. 


148 


long time 
solutions 


0.0020 
= Decreasing critical 
we heat flux 
© 
<= 
= 0.0015 
*) 
a) 
= 
— 
1S) 


0.0010 
---- Q",=25 kW/m’, T;, = 307 °C, 6=0.28 
—- a. Wp = 35 kWh’, T,, = 332 °C, >= 0.31 
9", = 50 kW/m’, T,, = 355 °C, 6 = 0.33 
Vr = 5 KW’, T,, = 377 °C, 6 = 0.35 


small time 


0.0005 solutions 


0.0000 
0 10 20 30 40 50 60 70 80 90 
Time, t {s] 


Figure 60 Effect of critical heat flux on char depth, dimensional plot. 


d) Heat of gasification, L 

Values for the heat of gasification of 100,000 J/kg, 1,500,000 J/kg, 2,000,000 J/kg and 
2,500,000 J/kg were selected in order to investigate its contribution to the burning rate 
and char depth (Figure 61 to Figure 65). Since the heat of gasification is the amount of 
energy required to vaporise the virgin wood, it would be expected that as it increases so 
the burning rate would decrease as shown in this analysis. However, the analysis also 
shows that as the heat of gasification is reduced, the small-time solution reaches a point 
where the burning rate begins at a higher value and immediately starts to decrease. Once _ 
this occurs, the small-time and long-time solutions for both the burning rate and char 


depth no longer intersect. 
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Figure 61. Effect of heat of gasification on burning rate, dimensionless plot. 
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Figure 62. Effect of heat of gasification on burning rate, dimensional plot. 
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Figure 63. Effect of heat of gasification on char depth, dimensionless plot. 


0.0030 


0.0020 


Char depth, 6, [m] 


0.0005 


0.0000 


if 
/ ee 


fA _-- longtime 
/ - solutions 


Decreasing heat 
of gasification 


~~ - L= 1,000,000 J/kg 
_ ~~. L= 1,500,000 J/kg 
__ L=2,000,000 J/kg 
__ L = 2,500,000 J/kg 


small time 
solutions 


0 10 20 30 40 50 60 70 80 90 
Time, t {s] 


Figure 64. Effect of heat of gasification on char depth, dimensional plot. 
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Figure 65. Effect of heat of gasification on char depth, dimensional plot with extended 
axes. . 


e) Apparent thermal inertia, koc 

In Figure 66 to Figure 69 the apparent thermal inertia was selected to be 

100,000 J?.m“Ks™, 200,000 J°.m“K”s", 300,000 J?.m*K”s" and 400,000 J?.m“K”s". 
In the analysis of the experimental data, the apparent thermal inertia was found to vary 
widely depending on the wood species and grain orientation. Since the specific heat 
capacity and thermal conductivity were calculated as a function of the thermal inertia 

(§ 2.7), they also varied as indicated. The analysis shows that as the thermal inertia 
increases, the model reaches a situation in which the small and long time burning rate 


solutions do not cross. 
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Figure 66. Effect of thermal inertia on burning rate, dimensionless plot. 
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Figure 67. Effect of thermal inertia on burning rate, dimensional plot. 
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Figure 68. Effect of thermal inertia on char depth, dimensionless plot. 
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Figure 69. Effect of thermal inertia on char depth, dimensional plot. 
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f) Thermal diffusivity, @ 

Although, as described in § 2.7, the thermal diffusivity is assumed to be constant, the 
dimensionless equations allow us to investigate its contribution to the model. Values of 
1.0x 10° m’/s, 2.1 x 10” m’/s, 3.2 x 10” m7/s and 4.3 x 10°’ m/s were graphed as shown 
in Figure 70 to Figure 73. Since the thermal diffusivity along with the thermal inertia are 
used to calculate the specific heat capacity and the thermal conductivity, the form of the 
curves are similar to those in which different thermal inertias are graphed as shown in 
part e) of this section. However, in this case, a decrease in the thermal diffusivity rather 
than an increase in the thermal inertia eventually leads to a situation where the small and 


long time burning rate solutions do not cross. 
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Figure 70. Effect of thermal diffusivity on burning rate, dimensionless plot. 
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Figure 71. Effect of thermal diffusivity on burning rate, dimensional plot. 
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Figure 72. Effect of thermal diffusivity on char depth, dimensionless plot. 
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Figure 73. Effect of thermal diffusivity on char depth, dimensional plot. 


g) Density, p 

The effect of density on the burning rate and char depth was investigated by choosing 
values of 400 kg/m’, 500 kg/m’, 600 kg/m? and 700 kg/m? (Figure 74 to Figure 77). The 
small-time burning rate solution is more sensitive to density and the long-time solution 
shows relatively little variation. In the work conducted by JFRO [48] it was found that 
the burning rate increased as the density decreased. The solutions to the integral model do 
not show this result but instead indicates that, although the peak burning rate is higher for 


lower densities, for long-times the burning rate is slightly higher for higher densities. 


In contrast to the burning rate, the variation of density has a greater effect on the long- 


time char depth solution with less dense materials charring quicker than more dense 
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materials. It is interesting to note that the transition between the small-time and long-time 


solutions is almost constant (Figure 77). 
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Figure 74. Effect of density on burning rate, dimensionless plot. 
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Figure 75. Effect of density on burning rate, dimensional plot. 
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Figure 76. Effect of density on char depth, dimensionless plot. 
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Figure 77. Effect of density on char depth, dimensional plot. 


h) Ignition time, tig 

Ignition times of 10 s, 20 s, 30 s and 40 s were selected to examine their contribution to 
the predicted burning rate and char depth (Figure 79 to Figure 81). In the dimensionless 
plots the curves fall on top of one another. The dimensional plot peat shows how the 
effect of increasing the ignition time shifts the burning rate curve along the x-axis but 


does not change the magnitude of the intersection point. 
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Figure 78. Effect of ignition time on burning rate, dimensionless plot. 
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Figure 79. Effect of ignition time on burning rate, dimensional plot. 
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Figure 80. Effect of ignition time on char depth, dimensionless plot. 
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Figure 81. Effect of ignition time on char depth, dimensional plot. 
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1) Char fraction , ¢ 

The dimensionless analysis (§ 4.1.6) shows that the char fraction only contributes to the 
long-time solution and this is shown in Figure 82. However, since the definition of the 
dimensionless burning rate M includes the char fraction, when the small-time 
dimensional burning rate is plotted, the char fraction makes a significant difference. As 


the char fraction approaches zero, the peak burning rate approaches that what would be 


given by a non-charring material 
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Figure 82. Effect of char fraction on burning rate, dimensional plot. 
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Figure 83. Effect of char fraction on burning rate, dimensionless plot. 
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Figure 84. Effect of char fraction on char depth, dimensionless plot. - 
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Figure 85. Effect of char fraction on char depth, dimensional plot. 


4.2 Experimental limitations and difficulties 
A number of limitations and difficulties were encountered during the course of the 


experimental programme. 


a) Scan rate 

The relatively slow scan rate used by the Cone Calorimeter may have resulted in aliasing 
problems particularly for obtaining the peak mass loss rate data. Unlike for a non- 
charring material, the peak mass loss rate of a charring material only occurs for a very 
short (almost instantaneous) time. If the peak mass loss occurred in-between scans then 
the actual value was not recorded rather a lower value may have been recorded (Figure 


86). 
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Figure 86. Error between actual and measured peak burning rate. 


In a few selected experiments the scan rate was increased to 1 or 2 scans per second to 
investigate whether this would enable a more accurate peak burning rate to be obtained 


although the results of this procedure were inconclusive. 


b) Thermocouples 

It was found that the use of thermocouples inserted into the sample could affect the mass 
loss rate measurements. The weight of the probes and the associated wires may have 
introduced an error and it was found that vibrations in the probes could result in very 
noisy mass loss data. For this reason, tests conducted later in the main ‘burning rate’ 
series were completed without any thermocouples inserted into the samples. As a result 
of the noisy mass loss rate measurements, the effective ate loss rate (or the effective 


burning rate) of the sample was obtained from Equation 1, using the measured rate of 
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heat release at each time interval and the average heat of combustion over the duration of 
the test. Since the instantaneous heat of combustion is obtained from the rate of heat 
release and the mass loss rate, the effect of the noisy mass loss rate data resulted also in 
noisy heat of combustion data and thus no average heat of combustion for a given test. In 
these cases the average heat of combustion calculated for all similar species of wood at a 


given orientation was used to obtain the effective mass loss rate. 


b) Back effects 

The integral model described in this study assumes that the sample is semi-infinite. As 
described in § 1.4.2, the samples necessarily had a finite thickness. For long duration 
exposures the thermal wave would effectively be ‘reflected’ at the back thus affecting the 
burning characteristics of the sample. Figure 83 shows an analysis for the predicted depth 
of the thermal penetration wave depth as a function of time using the long-time solution 


to the integral model (Equation 152) and the typical wood properties given in § 4.1.7. 


167 


Increasing incident 
heat flux, q”; 


Thermal penetration depth, 5 [m] 


Figure 87. Predicted depth of the thermal penetration wave using the long-term integral 
model solution. 


The integral model suggests that the thermal penetration wave reaches the back of the 
sample (i.e. a depth of 50 mm) in around 960 seconds (16 min). Thus, back effects may 
not be important up to this time. It should be note that the integral model developed by 


Moghtaderi et al. [31] successfully predicts the finite thickness effects. 


The finite thickness of the sample and the experimentally configuration also affected the 
thermocouple measurements. As the retainer frame was heated by the incident flux, it 
would transfer heat into the sides and back of the sample. This heat would penetrate into 
the sample and thus affect the internal temperature profiles. This effect was particularly 
noticeable with the thermocouple mounted at the back of the sample. The output of this 


thermocouple would often be higher than the thermocouple at 36 mm for the initial stages 
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of a test. Eventually, as the heat from the exposed surface penetrated through the sample, 


the 36 mm thermocouple would record temperatures above the back thermocouple. 


4.3 Data analysis and presentation 

A significant volume of data was recorded in this study for the burning rate of the 
samples and the accompanying temperature data. In addition the analysis of this data also 
generated a significant volume of information. In order not to present the reader with 
large numbers of plots in the main text, the complete set data are provided in the 
Appendix. For each sample tested, plots of the energy release rate; heat of combustion; 
temperature profiles; measured and predicted burning rate; derived and predicted thermal 


penetration depth and derived and predicted char depth are given where available. 


4.4 Cone Calorimeter data 
The energy release rate and instantaneous heat of combustion data for each sample tested 


are shown in the Appendix. 


4.4.1 Effective heat of combustion 
The average effective heat of combustion for an individual test was found as described in 


§ 1.4.2 and the overall average heat of combustion calculated for each species and grain 


orientation (Table 16). 


Table 16. Average effective heats of combustion for each species and grain orientation. 
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In comparison Tran & White [34] obtained values of between 10.5 kJ/g and 14.5 kJ/g for 


Redwood and between 9.0 kJ/g and 11.9 kJ/g for Red oak. 


4.4.2 Energy release 

In order to verify that the measurements recorded in this study were consistent, the 
average and peak energy release rates obtained from the oxygen calorimetry 
measurements were compared with data reported by Tran & White [34], Janssens [1] and 


Hopkins [30]. 


a) Douglas fir 
The peak energy release from the along grain samples compare well with J anssens [1] 
data (which were also tested in the same grain orientation) and the general trend in all 


three sets of data are similar. 
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Figure 88. Comparison of peak energy release rates from Douglas fir. 


b) Redwood 

Figure 89 shows that Janssens [1] obtained generally higher average energy release rates 
for Redwood compared to those obtained in this study. The peak energy release rates 
measured by Janssens fall in-between the along and across grain measurements recorded 


here (Figure 90). 
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Figure 89. Comparison of average energy release rates from Redwood. 
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Figure 90. Comparison of peak energy release rates from Redwood. 
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c) Red oak 

Figure 91 compares the average energy release from the Red oak samples with those 
quoted by Tran & White [34]. The along grain sample tested at 25 kW/m? is comparable 
with Tran & White’s data but further comparisons cannot be made since Tran & White 
did not measure the energy release rates at 75 kW/m’. 
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Figure 91. Comparison of average energy releases from Red oak. 


d) Maple 

Figure 92 shows the average energy release data for Maple. At an incident heat flux of 
75 kW/m’ the average values are almost identical for the two grain orientations. At 

25 kW/m’, there is a wider variation in the average energy release rate for the two tests 


conducted on the across grain samples. 
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Figure 92. Average energy releases from Maple. 


4.5 Char fraction 

4.5.1 Analysis 

In order to use the model to predict the burning rate of wood the char fraction needs to be 
calculated or measured. The char fraction for each sample tested in the experiments was 
found using Equation 5. Table 17 to Table 20 show the values used to obtain the char 


fractions for each specie. 
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Table 18. Sample dimensions used to obtain char fraction of Redwood. 
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0.0096 
0.0096 


Table 20. Sample dimensions used to obtain char fraction of Maple. 
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Figure 93. Char fraction against dimensionless irradiance. 


A plot of char fraction against dimensionless irradiance is shown in Figure 93. Assuming 


that the char fraction is related to the dimensionless irradiance such that 


g=a B 


where a and n are constants. By taking logs of both sides 


n 


1 
Ing = In| aj — 
: (5 
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By plotting Ing against nf 5) the values for a and n can be obtained from the gradient 


and intercept of the best-fit line. Figure 94 shows the In-In plot with a best-fit line. 
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Figure 94. In (char fraction) against In (dimensionless irradiance). 


Thus, the values for a and n were 0.74 and -0.64 respectively and the relationship 


between the char fraction and the dimensionless irradiance is 
1 —0.64 
g¢ = 0.74 (=) 
B 
Equation 156 


This relationship is plotted back on Figure 93 to show how it compares to the data. There 


is a fair degree of scatter of the data which may partly be due to the difficulties of 
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measuring the char depth. Clearly, this relationship may not be universal for all charring 


materials or even all species of wood but is applicable to this study. 


4.5.2 Measurements 

Three samples, all exposed for 75 minutes, were carefully examined and the actual char 
fraction obtained by measuring the mass of the char and the mass of the virgin wood. 
Figure 95 {a} shows sample 1ML5 prior to separating the char from the virgin wood. 
Note how it is not easy to determine the exact transition between the char layer and the 


virgin wood, which made the estimation of the volume difficult. 
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Figure 95. Half of samples of {a} 1ML5 and {b} 2DFX9i showing char layer and virgin 
wood. 
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0.000152 


0.000050 


0.000095 
0.000034 
0.0 
0.0204 
228 


0.000128 
0.000055 


1 337 606 576 
Le 


Table 21. Comparison of measured and calculated char fraction. 


Table 21 shows a comparison of the measured and calculated char fraction. Only with the 


10L2 sample is the comparison relatively close. It would have been desirable to have 


continued this comparison with additional samples, however, the samples were required 


by Schroeder [10] as part of his study and so were not available to the author. 


4.6 Burning rate 


4.6.1 Determination of flame heat flux and heat of gasification 


a) Burning rate analysis 


Equation 79 suggests that the net surface heat flux just after ignition is given by 


qf = 4, +qf-o(T! - TS) 


by neglecting the ambient temperature contribution which is small, this equation becomes 


ii =a +41-oF, 
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Equation 157 


Consider the case of a material with a steady state burning rate m/,,,,,. From Equation 


jf 3 


PIS Pécs( jp d 
2 ded 


“7}s + mn'c, (I, -T,)+ aed g’(T.)- eh) 


(ie 5, 


assuming that the surface temperature 7, of the burning material is at 7,, the thermal 


penetration depth 6, is 6, and the burning rate m” is mj,,,,, then 
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Equation 158 
Multiplying Equation 158 through by c,, 


Pans fart |= o, (t,)- eet) 


™(-¢) 6, 


Equation 159 
From the definition of ZL, Equation 159 can be written as 


2k, (L-AH,) 
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ee a 2k,L | 2k 
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substituting for 6, (§ 4.1.6) 
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| easy GT) _ G(T.) , a2(0,) 
(i-¢) AH, AH, L 


thus, recalling that 7, = 7;, the steady mass loss rate is given by 


» _(-9)9"(Z) 


™ steady ig 


Equation 160 
Combining Equation 157 with Equation 160 gives 
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Equation 161 


Using a linear regression line through the a plot of steady mass loss rate against incident 


heat flux, the slope can be expressed with the following equation 


fins Am" 
slope = Use wa 
Ad; 
Equation 162 
Thus 
re slope 
(1-4) 
Equation 163 


In addition, the value for the intercept of the best fit-line is where the incident heat flux 


q; equals zero. Therefore, Equation 161 can be rewritten as 


Tt eady = intercept = CG, -oT.’) 
Equation 164 
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By rearranging, we can calculate the flame heat flux gj, with 


., intercept - 
_ intercept-L on 


¥ (cae 


Equation 165 
assuming that the surface temperature 7; is at the ignition temperature Tig. 


The theory described above in Equation 164 and Equation 165 is similar to those 
developed for a non-charring material but with addition of the char fraction. For a 
charring material the mass loss rate does not reach a steady value at ignition but reaches a 
peak and then decays over time (§ 1.2.2). Thus, by plotting the peak mass loss rate 
against incident heat flux, an (over-)estimate of the heat of gasification and flame heat 


flux at the peak burning rate can be obtained from Equation 162 and Equation 165. 


Figure 96 to Figure 99 shows the graphs of the measured peak burning rate against 
incident heat flux for each species tested in the Cone Calorimeter. The heat of 
gasification and flame heat flux obtained from these in the along grain and across grain 


configurations are shown in Table 22 and Table 23. 


183 


30 


e Along grain 
© Across grain 
Along grain fit 
_ - - Across grain fit 


Peak mass loss rate, m"',,,, [g/s.m’] 


Incident heat flux, q"’, [kW/m] 


Figure 96. Peak burning rate as a function of incident heat flux for Douglas fir. 
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Figure 97. Peak burning rate as a function of incident heat flux for Redwood. 
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Figure 98. Peak burning rate as a function of incident heat flux for Red oak. 
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Figure 99. Peak burning rate as a function of incident heat flux for Maple. 
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Table 23. Calculated flame heat fluxes from peak burning rate data. 


b) Iterative approach 
An iterative approach, as suggested by Quintiere & Anderson [37], was used to adjust the 
flame heat flux and heat of gasification such that there was an overall best agreement 


between the model predictions and the experimental data. 


For each wood at a given orientation, the effective burning rate obtained from the Cone 
Calorimeter experiments were plotted in dimensionless terms and each individual 
dimensionless burning rate curve was plotted on a single set of axes (for e.g. Figure 96, 
similar results were obtained for the other woods). Since the definition of the 


dimensionless burning rate M is 


1 


(-9)a7(Z,) 


M =—— 
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initial values of the heat of gasification L and the flame heat flux gj, (used to find 47 (7. ) 


from Equation 81) were taken from the effective steady state burning rate analysis 


described above. 


Using the examination of the contribution of the flame heat flux and heat of gasification 
(§ 4.1.7) to the small-time and long-time solutions to the integral model as a guide, single 


values for these two variables were selected such that the following were achieved: 


12 The small-time solution gave a reasonable match with the initial mass loss rate at 
ignition, 
Il. The intersection of the small-time and long-time solutions was comparable with 


the peak experimental burning rate, 
iil. The long-time solution followed the decay portion of the experimental burning 


rate. 


It was found that adjusting the flame heat flux and heat of gasification to exactly match 
any one of the three criteria resulted in a poor comparison with the remaining two. Thus 
the selection of the flame heat flux and heat of gasification was based on obtaining a 


reasonable match with each of the three criteria. 


Clearly the values of the flame heat flux and heat of gasification obtained in this study are 
not intended to be exact for a particular species of wood at a given orientation. As already 


discussed in § 1.2.1, wood by its very nature varies from tree to tree and within any one 
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tree. In addition, the iterative process and the need to arbitrarily decide what constituted a 


good match means that the values given here are only representative. 
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Figure 100. Dimensionless burning rates of Douglas fir in the a/ong grain orientation 
using properties calculated in each individual test. 


0.20 


Table 24 and Table 25 show the final derived values for the heat of gasification and flame 


heat flux for each wood in the two grain configurations. 


Table 24. Final derived heats of gasification from iterative analysis. 
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Table 25. Final derived flame heat fluxes from iterative analysis. 


In the study by Rhodes & Quintiere [46] and also by Hopkins [30] it was found that the 
flame heat flux for burning PMMA was approximately constant at 37 kW/m?. Hopkins 
also measured the flame heat flux as a function of time for Redwood and Red oak. For 
Redwood, peak flame heat fluxes between 32 kW/m” and 53 kW/m’ with an average of 
44 kW/m’ were recorded by Hopkins compared to an average of 34 kW/m’ obtained in 
this study. For Red oak, Hopkins recorded values between 32 kW/m? and 52 kW/m’ with 
an average of 42 kW/m’ compared to an average of again 34 kW/m’ obtained in this 
study. In comparison, for thermoplastics, Hopkins obtained values between 14 kW/m’ 
and 37 kW/m’. The peak flame heat fluxes obtained by Hopkins are on average higher 
than those obtained using the iterative approach. However, it should be noted that in 
many of the experiments conducted by Hopkins, the flame heat flux would reduce over 
time presumably because, as observed in § 1.4.2, the flames became smaller as the char 
layer formed. We might conclude that the flame heat fluxes derived here are more akin to 
the average flame heat flux over the duration of the test and thus would be expected to be 


less than the peak value. 


Figure 101 shows the dimensionless burning rate for Douglas fir, along grain, for the four 
incident heat fluxes using the final heat of gasification and flame heat flux given in Table 


24 and Table 25. The chart is a combination of the individual burning rates for each 
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sample, averaged and smoothed so as to obtain the general form of the experimental data. 
Comparison with Figure 49 shows how the experimental data compares with the integral 
model solutions. The peak dimensionless burning rate at low heat fluxes is greater than at 
high heat fluxes in both cases. The relative offset in the peaks as a function of time is 
more difficult to identify in the experimental data. At long-times the burning rates at the 
various heat fluxes run parallel although there is some minor fluctuation in the 


experimental data. Similar results were obtained for the other species. 
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Figure 101. Combined dimensionless burning rates of Douglas fire in the along grain 
orientation using derived properties. 


Figure 102 shows the dimensionless char depths for Douglas fir, along grain. Comparison 
with Figure 51 shows similarities between the experimental data and integral model 


solutions. In general, the char depths run parallel leading to an inflection point at around 
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T= 0.1 where the gradient of the lines decrease but continue to run parallel. The relative 


order of the experimental data does not compare as well with the integral model 


solutions. For example, Test 1DF3 (25 kW/m”) shows a higher dimensionless char depth 


than the other tests whereas the integral model predicts lower dimensionless char depths 


for lower incident heat fluxes. 
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Figure 102. Dimensionless char depths of Douglas fire in the a/ong grain orientation 


using derived properties. 


4.6.2 Comparisons with model predictions 


The Appendix shows the comparison between the measured effective mass loss rate (as 


discussed in § 4.2) and the small-time and long-time integral model solutions. To obtain 


the solutions, the analysis used the following properties and relationships: 
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1) the heat of gasification and flame heat flux derived from the iterative approach 
(Table 24 and Table 25), 

il) the critical heat flux (Table 11), ignition temperature (Table 12) and thermal 
inertia (Table 13) obtained from the ignition analysis, 

lil) the average density (Table 3) and thermal diffusivity (Table 4), 

iv) the convective heat transfer relationship (Equation 67) and 


v) the char fraction relationship (Equation 156) 


For each sample, the burning rate is shown in both dimensionless and dimensional forms 
with common scales for the x-axis and y-axis to allow for easy comparison. Examination 


of the burning rate curves found the following points of note: 


a) Douglas fir 

The peak burning rate obtained by the intersection of the small-time and long-time 
solutions is greater than the measured data at high incident heat fluxes. For the along 
grain samples the long-time burning rate is lower than the measured data however a 
better match is achieved for the across grain samples. ne discussed in § 4.6.1, by 
changing the properties we could get a better match between the long-time solution and 
the experimental data but at the expense of the match between the peak burning rates. 
Thus the results shown for this and the other woods is a compromise between these 
conflicting requirements. The results from Test 1DFX3 are discussed in more detail in 


§ 4.7. 
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b) Redwood 
For reasons not determined, the energy release data from Test 1RL1 shows an offset. For 
along grain samples the long-time burning rate is lower than the measured data but, as for 


Douglas fir, the across grain results compare well. 


c) Red oak 
Both grain orientations exhibit a reasonable general match between the measured burning 


rate and the integral model. 


d) Maple 

The across grain samples show a good match between the measured and predicted 
burning rate. The along grain samples again have a lower burning rate from the integral 
model compared with the measured data. Only limited data was available for Test IML5 
as the Cone Calorimeter’s data acquisition system crashed 3:50 after the start of the test. 
The system was restarted at 5:45 by which point the sample had just started to burn. 
Processing the data files for this test was problematic and thus some of the analysis has 


not been completed for this sample. 


In several cases, it appeared that the long-time solution for the burning rate using the 
derived properties gave an under-estimate compared to the measured data. As shown in 
§ 4.1.7, by decreasing the heat of gasification, the long-time burning rate can be shifted 
upwards. However if the heat of gasification is decreased sufficiently, it is found that the 


small-time and long-time solutions no longer intersect thus providing a solution that 
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departs from what is actually measured. Again, this demonstrates the difficulty in 


determining the properties whilst maintaining an overall reasonable match with reality. 


4.6.3 Peak burning rate 

As discussed in § 4.1.7, the small-time solution of the integral model eventually reaches a 
steady-state level equivalent to what would be expected from a non-charring material 
(Figure 42). The peak burning rate from a charring material, whether measured or 
obtained from the intersection of the small-time and long-time solutions, is always less 


than the small-time steady-state level as shown in § 4.1.7. 


By calculating the small-time steady-state value and obtaining the peak burning rate 
either from experiments or the integral model, the ratio 


Ta eat 


Ti veady 
Equation 166 

can be found for each sample at its incident heat flux exposure condition. Figure 103 and 
Figure 104 show the average ratios for each wood species at the given orientations and 
incident heat fluxes. The general form of the curves is similar with the ratio higher at the 
low incident heat fluxes and decreasing as the heat flux increases. In general, the curves 
obtained from the integral model compare well. For across grain Maple the ratio remains 
almost constant and for along grain Douglas fir and Redwood the ratio increases at 35 


kW/m” compared to all species apart from the across grain Maple. 
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Figure 103. Ratio of charring peak burning rate with non-charring steady-state burning 
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Figure 104. Ratio of charring peak burning rate with non-charring steady-state burning 
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An overall set of average values can be obtained from the results of the above analysis of 


the charring peak burning rate with non-charring peak burning rate (Figure 105). In this 


case, where no ratio for the experimental data exists for a particular wood (e.g. Maple, 


across grain at 35 kW/m? etc.) then the ratio is found by interpolation. The chart indicates 


the reasonable comparison between the integral model and the experimental data over the 


range of incident heat fluxes. 
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Figure 105. Average ratios of charring peak burning rate with non-charnng peak burning 


rate from integral model solutions and experimental data. 


Thus, an empirical relationship can be obtained for estimating the peak burning rate from 


the integral model solution. From Equation 160, the steady burning rate is given by 


.» _(-9)4i(Z,) 


N steady L 
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Using a cubic polynomial fit to the integral model ratios as shown in Figure 105, the peak 


burning rate factor /7 is given by 
TI = -0.000004[g’]’ + 0.00083[g"] - 0.06g"+1.96 


Equation 167 
thus, from Equation 166, the peak burning rate is simply 


" 


. ° n 
MN peak ves TT ready 


Equation 168 


Figure 106 shows a comparison between the calculated peak burning rate using Equation 


167 and Equation 168 with the experimentally measured peak burning rates. The values 
for m‘,..4, Were obtained using the values for the heat of gasification and flame heat flux 


given in § 4.6.1 (Table 24 and Table 25). The char fraction ¢ was found from Equation 


156. 
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Figure 106. Comparison of measured peak burning rate and calculated burning rate using 
the ratio of peak to steady burning rate relationship. 


The empirical relationship for the prediction of the peak burning rate shows that, 
although the trend is consistent, for higher measured values the relationship under- 


estimates the peak burning rate. 


4.7 Burning rate of char 

During the tests conducted in this study it was observed that flaming would eventually 
cease and instead the char at the surface of the material would be oxidised away by the 
incident heat flux. The best example of this was in Test 1DFX3 (Douglas fir, 75 minutes 
at 75 kW/m’) where the rate of heat release and heat of combustion data clearly shows 


the transition from wood burning to char oxidation (Figure 107). 
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Figure 107. Rate of heat release and heat of combustion for Test 1DFX3 showing back 
effect and char oxidation phases. 


The heat of combustion data exhibits a higher degree of fluctuation during the char 
oxidation phase due to the relatively low instantaneous mass loss rates. From Equation 1 
it can be seen that small values for the instantaneous mass loss have a significant effect 
on the calculated instantaneous heat of combustion. By taking an average of the heat of 
combustion during the char oxidation stage of the test and neglecting the high and low 
peaks, the heat of combustion of the char was found to be 35.5 MJ/kg. This value 
compares well with the average heat of combustion for char of 34.3 MJ/kg quoted by 
Drysdale [2]. By using the average heat of combustion for char and the measured rates of | 


heat release the average ‘burning’ rate of the char was found to be 0.0015 kg/s.m?. 
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Figure 107 also clearly shows the back effect (§ 4.2) where the burning rate of the sample 


increases as the thermal wave is reflected at the rear of the sample. It is interesting to note 


that the increase in burning rate appears to begin at around 20 minutes. This result is 


comparable to the value of around 16 minutes for the thermal wave to reach a depth of 


50 mm found in the analysis of the long-time solution for the integral model (§ 4.2). 


4.8 Charring rate and char depth 


4.8.1 Background 


In this section the char depth is analysed using the integral model solutions, the 


thermocouple measurements, the mass loss data and the measured char depth at the end 


of the test. These aspects are discussed and compared below. 


4.8.2 Mass loss rate data 


From the derivation of governing equations, Equation 125 states that 


m’ do, 


(=) ae 
thus rearranging 


do, sm" 
dt p,(1-) 


and integrating over time t we obtain 


on 


m 


0.= 
*  p(l-¢g) 


t 
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Equation 169 


Equation 170 


Since the experimental data provides the char depth at the end of each test, we can 


compare the measured char depth and the calculated char depth from Equation 170 using 


the average mass loss rate measured in the Cone Calorimeter and ¢ equal to the duration 


of the test ty. Figure 108 shows a plot of calculated char depth against measured char 


depth. At earlier times (25 minutes), the calculated and measured data match fairly well 


though the mass loss rate relationship appears to over-predict slightly. However, at later 


times (75 minutes) the integral model fails to successfully predict the char depth 


compared with the measured data. It is likely that this is due to a back face effect where 


in the experiments heat is lost through the back face and/or the sample is completely 


charred whereas the integral model treats the sample as having an infinite depth that can 


continue to char indefinitely. 
0.10 
| @ Douglas fir - along gran 
0.09 © Douglas fir - across grain 
a Redwood - along grain 
0.08 4 pb Redwood - across grain 
a Red oak - along gram 
= 007 | a Red oak - across grain 
Pe @ Maple - along grain 
= 0.06 © Maple - across grain 
2 
& 0.05 
s 
= t= 25 mins oO 
= e 
Se 0.03 
0.02 i ° 
ieee 
0.01 3 #0 


oe 
<-> me 
on Te. 


oo 


t;= 75 mins _ 
—> » 
‘ 
cy 
oe 
at 
eee 
‘ ae 
eee oe \ a Equality lme | 
oe atu a ' 
--~ ba ee H 
a faz - 
Bete a ‘te oe i 
8 o ig Hi 
ed ‘ 
e eg (ahs - ed ae ie , 
‘. 
e@ ~2 Q. 3 ° = at / ‘ a 
- Z . 
ae i | Ke “A 
f ‘s, : 
eo wi as wa 
me va Sec re) Bosc 
a th ~ -" 


=. SS 
= 
2 SS eeeceno="~ 


Measured char depth, 5, [m] 


Figure 108. Comparison of measured and calculated char depth. 
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It is interesting to note that Tran & White [34] suggest that the charring rate can be found 


from 


Equation 171 


assuming that the mass of the char is negligible and ignoring the mass loss in advance of . 
the char front. Equation 171 is equivalent to the charring rate equation, Equation 169, 


obtained from the integral model with the difference being that the integral model 
includes the effect of the char fraction. In § 4.5.1 we show that as gj >>q? the ¢ > 0, 


thus Equation 169 and Equation 171 would give the same result. 


4.8.3 Final char depth from integral model solutions 
A similar analysis as that shown in § 4.8.2 can be performed using the predicted char 
depth from the long-time solution given by the integral model. The dimensionless char 


at 


depth is given by Equation 154. From t= 5 the equivalent dimensionless time 7 for 


25 minutes can be obtained and hence the char depth calculated. Similarly the measured 
char depth and the char depth from the mass loss data can be made dimensionless since 


O4 


AC 
¢ 6, 


Thus a comparison can be made by plotting the dimensionless char depths from the 
integral model, the measured data and the mass loss data against each other (Figure 109 


to Figure 112). In each figure, three plots are shown; the measured depths versus the 
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integral model solution (line 1, using + symbols), the mass loss versus integral model 


solution (line 2, using - symbols) and mass loss versus measured depths (line 3, using X 


symbols). 
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Figure 109. Comparison of char depths at 25 minutes for Douglas fir. 
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Figure 110. Comparison of char depths at 25 minutes for Redwood. 
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Figure 111. Comparison of char depths at 25 minutes for Red oak. 
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Figure 112. Comparison of char depths at 25 minutes for Maple. 


The analysis shows, as found in § 4.8.2, the mass loss and measured depths compare 
well. The data (line 3) lie close to the equality line and follow the equality line over the 


range of depths obtained. 


The comparison between the integral model and the measured data (line 1) shows that at 
low dimensionless char depths (between 1.0 and 2.0) the two compare well but as the 


depth increases the integral model predicts shallower depths than measured. 


4.8.4 Char depth using thermocouple data 
In addition to using the thermocouple data to obtain the thermal penetration wave 
(§ 4.9.1), the measurements were used to map the travel of the pyrolysis front (i.e. char 


depth) as it moved through the samples. 
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In § 4.1.1 it is assumed that the vaporisation temperature of the volatiles 7, in the 
pyrolysis zone is equivalent to the ignition temperature Tig. Thus, if the time at which the 
thermocouples first record a temperature at the calculated ignition temperature for a 
species of wood at a given orientation (Table 12) then the char depth can be obtained. In 
the studies by Schaffer [17] it was found that the transition temperature for wood to 
become char is around 288 °C regardless of species. Therefore, in this study, the char 
depth using the thermocouple measurements has also been obtained by using a critical 
temperature of Tar = 288 °C in addition to the calculated ignition temperatures. Figure 


113 shows an example of the method used to obtain the char depths. 
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Figure 113. Determination of char depth using wood ignition temperature or char 
ignition temperature (Test 1RL2). 
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The char depth can be examined in terms of dimensionless variables for each species of 
wood at given incident heat fluxes. The dimensionless depth of each thermocouple probe 


can be obtained from 


Equation 172 
where 6,- is the depth of the thermocouple and the time at which it reaches the ignition 


temperature ¢(Tjz) from 


Equation 173 
The integral model gives Equation 143 for the small-time char depth, 


A, * 2M, (r-z,) 


and Equation 154 for the long-time char depth, 


A, oy {22 Joe ~1)(r-7,) 


Plots for the char depth using the thermocouple data are compared with the theoretical 
predictions from the integral model are presented below. In each plot a trend line is 
drawn through the thermocouple data as a means to compare the experimental data with 


the theoretical curves. 


a) Douglas fir 
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Figure 114. Dimensionless char depth using thermocouple measurements for Douglas fir 
at 25 kW/m’. 
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Figure 115. Dimensionless char depth using thermocouple measurements for Douglas fir 
at 50 kW/m’. 
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Figure 116. Dimensionless char depth using thermocouple measurements for Douglas fir 
at 75 kW/m’. 


Figure 114 to Figure 116 show the dimensionless char depth plots for Douglas fir. The 
trend for the thermocouple data at 25 kW/m? differs from the theory simply because only 
two depths were obtained from the temperature data. The plots for 50 kW/m? and 

75 kW/m? thermocouple data show reasonable consistency between tests particularly for 
the across grain configuration. Comparison of the thermocouple data with the theoretical 


predictions show a fair qualitative match. 


b) Redwood 
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Figure 117. Dimensionless char depth using thermocouple measurements for Redwood 
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Figure 118. Dimensionless char depth using thermocouple measurements for Redwood 


at 50 kW/m’. 
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Figure 119. Dimensionless char depth using thermocouple measurements for Redwood 
at 75 kW/m’. 


Figure 117 to Figure 119 show the dimensionless char depth plots for Redwood. The 
plots for the thermocouple data show good consistency in both grain configuration. 
Comparison of the thermocouple data with the theoretical predictions show a good 
qualitative match compared with the Douglas fir. In all three cases the dimensionless char 
depth for the across grain configuration is greater than the along grain configuration for 
both the thermocouple data and the integral model predictions. The general trend for the 
thermocouple data follows the model predictions albeit the thermocouple data gives 


generally higher dimensionless char depths compared with the theoretical values. 


Similarly, dimensionless plots could be made for Red oak and Maple and also plots using 


the critical char temperature of 288 °C. 
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4.8.5 Overall results 
The char depth obtained from the mass loss data, the integral model solutions and the 
thermocouple data for each wood sample tested is shown in the Appendix. Concluding 


observations regarding the results for each wood species are given here: 


a) Douglas fir 

For the along grain samples, with the exception of Test 1DFL3 (25 kW/m? for 

25 minutes), the integral model follows the mass loss data up to around 8-10 minutes at 
which point the two curves intersect. Thereafter, the integral model predicts shallower 
char depths than given by the mass loss data. This result agrees with the findings 
discussed in § 4.8.3. For Test 1DFL3 the short-time solution prevails for almost the 
whole test. In general, the thermocouple data compares reasonable well with the mass 
loss and integral solutions with the 7;, analysis comparing better in some cases and the 


T char analysis 1n others. 


For the across grain samples the integral model and mass loss data show similar trends 
with the long-time solution and mass loss data almost coinciding at the termination of 
each test. The 7;, thermocouple measurement analysis matches the integral model 
reasonably well. In Test 1DFX3 the thermocouple data follows the mass loss data 


remarkably well (Figure 187) as the end effects become significant. 


b) Redwood 
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For the along grain samples the long-time solution and mass loss data only compare well 
for the first few minutes and thereafter the integral model predicts shallower char depths. 
With the exception of Test 1RL1, the 7;, thermocouple data closely follows the mass loss 
data. In Test 1RL1 (25 minutes at 75 kW/m”), there appeared to an offset with the 


measured mass loss data that affect the analysis. 


The integral model and char depth from mass loss rate measurements compare well for 
across grain samples. The thermocouple data follows the same trends with the Tena, 


analysis performing better than the Jj, analysis. 


c) Red oak 

There is a good match between the integral model and the mass loss analysis of the char 
depth. The match between the char depth and the thermocouple data is fair. Similar to 
Test 1DFX3, in Test 1OL2 the thermocouple measurements again follow the mass loss 
rate at long times. There was only limited data for Test 1OL3 as the 12 mm and 24 mm 


deep thermocouples failed early in the test. 


d) Maple 

The along grain integral model predictions and mass loss rate char depths compare well 
however the thermocouple data gives deeper char depths in comparison. For across grain, . 
the T-nar temperatures match the integral model for an incident heat flux of 25 kW/m? 


whereas the Ti, data gives deeper char depths. For an incident heat flux of 75 kW/m’ the 
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integral model is comparable to the mass loss data but the thermocouple data does not 


compare well. 


4.9 Temperature measurements 
4.9.1 Thermal penetration depth 
From the integral model we can determine the depth of the thera penetration wave 
prior to ignition from Equation 27 
O’= 24/28, 
pe|1-B 
Once ignition has occurred, the thermal penetration wave follows Equation 151 from the 


small-time solution, 


saa, +6 iH 
AH 


v 


Je —A, Ne _ re 
and Equation 152 from the long-time solution, 


2 
Ax NAG +12(r-7, 


The measurements from the thermocouples inserted into the samples were used to map 
the thermal penetration wave and Figure 120 shows the analysis method. A temperature 
difference AT between ambient and a given temperature rise was selected. The 
thermocouple measurements at each depth were examined and the first time at which the 
measured value reached or exceeded AT was noted. Thus the time for the thermal wave to 
penetrate through the depth of the sample could be obtained as a function of the depth of 


each thermocouple. 
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Figure 120. Detail of temperature measurements for Test 1DFX4. 


It must be recalled that, as described in § 4.2, at greater depths and longer times the 
temperature measurements are influenced by the heat penetrating through the rear of the 
sample from the hot sample retainer. Figure 120 also shows how the 50 mm (back) 
thermocouple would record higher temperatures during the initial stages of a test due to 
the hot retainer frame. As the test progresses, the temperature from the 36 mm 


thermocouple continues to rise and eventually overtakes the back thermocouple. 


In theory, the moment the top surface of the sample is exposed to the incident heat flux, 
the rise in temperature is transmitted through the complete depth of the material. In 
practice this is not the case and even if this were true it would not be possible to measure 
the minute temperature changes at the locations remote from the exposed surface. Thus a 


temperature rise was selected that was relatively small so that the progress of the thermal 
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penetration wave through the sample was measured, but not so small that the inherent 


measurement error from the thermocouples and other factors became important. 


Figure 121 shows the measured thermal penetration depth using four temperature rises of 
2 °C, 5 °C, 10 °C and 15 °C and the theoretical values for Test 1DFX4 (25 kW/m’). It 
clearly shows how the results for the back thermocouple give erroneous results for the 
time for the thermal penetration wave to reach a depth of 50 mm. For a AT of 2 °C the 
time is only 42 s which is considerably sooner than the time of 170 s for the thermal wave 
to reach 36 mm. As AT is increased, the back thermocouple begins to give results that are 
more appropriate and by a AT of 15 °C the result is closer to what would be expected. 
Because of this problem with the back thermocouple, it was normally neglected from the 


analysis of the thermal penetration wave. 
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Figure 121. Determination of thermal penetration depth using temperature rises (AT) 
measured by the thermocouples (1DFX4). 


Figure 122 shows a similar analysis for Test IDFL1 (75 kW/m’). At small times, i.e. at 
the top-most thermocouple, there is little difference between the four temperature rises 
and the theory. As the thermal wave penetrates the sample, using the different 
temperature differences begins to have a greater effect on the match between the 
measurements and the theory. At 24 mm the 2 °C temperature difference shows the best 
match with the theory whereas at 36 mm the 15 °C temperaturé'différence shows the best 


match with the theory. 
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Figure 122. Determination of thermal penetration depth using temperature rises (AT) 
measured by the thermocouples (1 DFL1). 


Furthermore, the theory does not account for the effects of moisture in the sample. Which 
is evident in Figure 123. Some heat is required to vaporise the moisture which results in 
the plateaus that can be seen in the temperature measurements at 100 °C. This process 


effectively slows the thermal penetration wave as it moves through the sample. 
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Figure 123. Temperature profiles for Test 1DFL1 showing plateaus at 100 °C. 


As a result of this analysis and in order to obtain an appropriate approximation of the 
thermal penetration wave a temperature difference of 5 °C was selected for the 


comparison between the measured and theoretical temperature profiles. 


Plots for the thermal penetration depths in each sample in which thermocouples were 
utilised are shown in the Appendix. For all of the tests in which thermocouples were 
used the match between the measured data and integral model solutions is excellent to a 
depth of 36 mm. Only for the across grain Redwood samples do the thermocouples and 


integral model not compare so well. 
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4.9.2 Surface temperature 


Equation 153 gives the temperature of the exposed surface after a ‘long’ time where 


g" sys 1/4 
Oo 


If it is assumed that the maximum temperature recorded by the top-most thermocouple is 
equivetent to the surface temperature (i.e. the wood surface and the wood or char at a 
depth of 4 mm are at equilibrium) then the measurements and Equation 153 can be 
compared. As discussed in § 1.4.2 and 4.6.1, the flame on the buming surface of the 
sample would not remain constant but would gradually reduce and eventually extinguish. 


Thus Equation 153 gives an estimate of the maximum long time surface temperature. If 
it assumed that the flame heat flux at long time is negligible i.e. gj, ~ 0, then Equation 


153 becomes 


Equation 174 


and thus this is an estimate of the minimum long time surface temperature. The 
maximum temperature was obtained from each test in which thermocouples were utilised 
and plotted against incident heat flux (Figure 124 and Figure 125). The figures also show 
the theoretical maximum and minimum long time surface temperatures given by 


Equation 153 and Equation 174 using the flame heat fluxes derived in § 4.6.1 (Table 25). 
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Figure 124. Measured and theoretical long time surface temperatures for Douglas fir. 
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Figure 125. Measured and theoretical long time surface temperatures for Redwood. 


221 


Temperature, T, [°C] 
nN 
S 
4 


300 | ——— Theoretical without flame 
| — Theoretical with flame (along grain) 
200 4 w Thermocouple data (along grain) 


_ — ~ Theoretical with flame (across grain) 
1). Thermocouple data (across grain) 


0 10 20 30 40 50 60 70 80 
Incident heat flux, q", [kW/m’] 


Figure 126. Measured and theoretical long time surface temperatures for Red oak. 
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Figure 127. Measured and theoretical long time surface temperatures for Maple. 
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For Douglas fir and Redwood it can be seen that the thermocouple measurements 
generally fall in-between the two theoretical curves. For the Douglas fir along grain 
surface temperatures at high heat fluxes the thermocouple measurements fall below the 
theoretical minimum temperature curve. Similarly for Red oak and Maple, the 
thermocouple measurements generally fall in-between the two theoretical curves at 


25 kW/m” but lie below the theoretical curves at 75 kW/m’. 
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5. 


FURTHER WORK 


This study presents a relatively simple model describing the ignition and burning 


characteristics of charring materials and methods to obtain the properties required to use 


the model. Suggested further work is given here although some of these may have only a 


minor effect on the theoretical predictions. 


a) 


b) 


c) 


The current integral model assumes that the char layer remains in place throughout 
the burning duration. As shown in this study, the exposed surface of the char layer 
recedes as it is oxidised by the incident heat flux. Thus the model could be modified 
to account for the oxidation in terms of an additional mass loss rate from the sample 
and thus a reduction in the overall height of the sample 

The integral model does not explicitly include the effect of any free moisture in the 
sample. As discussed in § 2.3, the moisture content of wood will influence the 
ignition properties. However, this moisture effect is effectively accounted for in the 
apparent thermal inertia obtained from the time to ignition data. Moisture also plays a 
part in the burning rate of wood and, as discussed in § 4.9.1, the thermal penetration 
into the material. Again the effect of the moisture is, at least partly, accounted for by 
the derivation of the effective heat of gasification of the material. 

The processes involved in the ignition of wood at lower incident heat fluxes clearly 
warrants further investigation. The results from this study have shown that the 
mechanism for ignition appears to change from a ‘bulk’ surface ignition to a localised 
glowing ignition. The ignition temperature of wood is almost constant at high 
incident heat fluxes but appears to initially rise as the flux 1s reduced before falling to 


considerably lower values as the glowing ignition process becomes predominant. 


224 


d) 


Finally, there appears to be an almost constant transitional incident flux at which the 
ignition of wood as a function of grain exhibits a cross-over between the relative 
times to ignition. 

As observed in § 1.4.2 and measured by Hopkins [30] (§ 4.6.1), the flame from a 
burning charring material such as wood does not remain constant as assumed in the 
integral model. The size of the flame, its heat flux back to the surface of the material 
and thus the mass flow rate of the volatiles are all interdependent. We might therefore 


want to consider modification of the model to account for this process. 


CONCLUSIONS 
The integral model for the time to ignition gives good agreement with experimental 
data at high incident heat fluxes (greater than ~20 kW/m?) 
A low estimate of the critical heat flux for piloted ignition can be obtained from the 
time to ignition data using the intercept along the x-axis of a linear extrapolation of a 
plot of 1/ aan against incident heat flux. Comparison with the integral model theory 
shows that this intercept value needs to be modified by a constant factor to obtain a 
better estimate of the critical heat flux. 
An estimate for the ignition temperature of wood can be obtained from the critical 
heat flux derived from the ignition time measurements. 


The apparent thermal inertia of a material can be obtained from the slope of a linear 


extrapolation of a plot of 1/ Vig against incident heat flux. 


eRe 


The ignition and burning rate of wood depends on many factors including the species, 
grain orientation, moisture content, exposure conditions and the inherent variability of 
wood as a natural material. 

The integral model for the burning rate of a charring material combines a small-time 
and a long-time solution. It may be possible to use the transition between the two 
solutions to obtain an estimate of the peak burning rate of the material although the 
results obtained do not compare particularly well with the measurements. 

The effective heat of gasification and flame heat flux can be obtained from burning 
rate experiments conducted in the Cone Calorimeter. An iterative approach is 
necessary to obtain appropriate values of the heat of gasification and flame heat flux 
which best match the measured data. The choice of the properties is a balance 
between conflicting requirements so as to match the small-time solution, the peak 
burning rate and the long-time solution and is therefore based on a certain degree of 
interpretation. 

The char fraction of wood has been found to be a function of the ratio of the incident 
heat flux and the critical heat flux. The function may not be universal for all charring 
materials or even all species of wood. 

Once suitable properties have been derived, the integral model’s solutions for the 
burning rate qualitatively compares well with the measured data. The model predicts 
both the initial growth and the subsequent decay. 

Thermocouples embedded into a burning sample can be used to estimate the thermal 
penetration depth, char depth and the surface temperature after a long time. 


Comparison between the thermocouple data and the integral model shows reasonable 
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agreement particularly for the thermal penetration depth and the long-time surface 
temperatures. Prediction of the char depth using the thermocouple data is not 
particularly reliable. 

Prediction of the char depth can be made using the average buming rate. The 
predictions of char depth based on the measured mass loss rate and the measured char 
depth after an exposure of 25 minutes compare well. For an exposure of 75 minutes, 
the match : poor and this may be due to the finite thickness of the samples. 

The mechanism for the ignition of wood at low heat fluxes close to the critical heat 
flux appears to be different from that at high heat fluxes. At low heat fluxes, localised 
smouldering of the wood may increase the energy input at that point and thus lead to 


a localised ignition. 


“Science is not the easy road to simple answers, but a torturous path to more complex 
questions.” 


“And this thumping in our chest said ... there is no rest. 
We would hang together 
or they’d hang us one by one” 


A Fire is Burning, Oysterband (Holy Bandits, 1993) 
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Test (IDFL1), Douglas Fir, along grain at 75 kW/m’ for 25 minutes. 
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Figure 128. Rate of heat release and heat of combustion (1DFL]1). 
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Figure 129. Temperatures measured in sample (1DFL1). 
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Figure 130. Comparison of dimensionless burning rate using derived properties for 
species and orientation (1DFL1). 
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Figure 131. Comparison of burning rate using derived properties for species and 
orientation (1DFL1). 
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Figure 132. Comparison of measured and calculated thermal penetration depth (1DFL1). 
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Figure 133. Comparison of estimated char depth from experimental data and the 


calculated char depth (1DFL1). 
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Test (1DFL2), Douglas Fir, along grain at 75 kW/m? for 25 minutes. 
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Figure 135. Temperatures measured in sample (1DFL2). 
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Figure 136. Comparison of dimensionless burning rate using derived properties for 
species and orientation (1 DFL2). 
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Figure 137. Comparison of burning rate using derived properties for species and 
orientation (1DFL2). 
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Figure 138. Comparison of measured and calculated thermal penetration depth (1DFL2). 
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Figure 139. Comparison of estimated char depth from experimental data and the 


calculated char depth (1DFL2). 
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Test (1IDFL3), Douglas Fir, along grain at 25 kW/m” for 25 minutes. 
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Figure 140. Rate of heat release and heat of combustion (1DFL3). 
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Figure 141. Temperatures measured in sample (1DFL3). 
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Figure 142. Comparison of dimensionless burning rate using derived properties for 
species and orientation (1DFL3). 
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Figure 143. Comparison of burning rate using derived properties for species and 
orientation (1DFL3). 
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Figure 144. Comparison of measured and calculated thermal penetration depth (1DFL3). 
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Figure 145. Comparison of estimated char depth from experimental data and the 


calculated char depth (1DFL3). 
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Test (1DFL4), Douglas Fir, along grain at 75 kW/m’ for 25 minutes. 
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Figure 147. Temperatures measured in sample (1 DFL4). 
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Figure 148. Comparison of dimensionless burning rate using derived properties for 
species and orientation (1DFL4). 
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Figure 149. Comparison of burning rate using derived properties for species and 
orientation (1DFL4). 
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Figure 150. Comparison of measured and calculated thermal penetration depth (1DFL4). 
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Figure 151. Comparison of estimated char depth from experimental data and the 
calculated char depth (1DFL4). 
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Test (1DFL6), Douglas Fir, along grain at 50 kW/m? for 25 minutes. 
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Figure 152. Rate of heat release (1DFL6). 
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Figure 153. Temperatures measured in sample (1DFL6). 


241 


——— Expermental data (1) 
~ —— Small-time solution (2) 
— —Long-time solution (3) 


Dimensionless burning rate, M [-] 


0.00 0.10 0.20 0.30 0.40 0.50 0.60 0.70 0.80 0.90 1.00 
Dimensionless time, t - t,, [-] 


Figure 154. Comparison of dimensionless burning rate using derived properties for 
species and orientation (1DFL6). 
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Figure 155. Comparison of burning rate using derived properties for species and 
orientation (1DFL6). 
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Figure 157. Comparison of estimated char depth from experimental] data and the 


calculated char depth (1 DFL6). 
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Test (1DFL7), Douglas Fir, along grain at 50 kW/m’ for 25 minutes. 
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Figure 158. Rate of heat release and heat of combustion (1 DFL7). 
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Figure 159. Temperatures measured in sample (1DFL7). 
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Figure 160. Comparison of dimensionless burning rate using derived properties for 
species and orientation (1DFL7). 
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Figure 161. Comparison of burning rate using derived properties for species and 
orientation (1DFL7). 
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Figure 162. Comparison of measured and calculated thermal penetration depth (1DFL7). 
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Figure 163. Comparison of estimated char depth from experimental data and the 


calculated char depth (1DFL7). 
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Test (1DFL8), Douglas Fir, along grain at 50 kW/m’ for 25 minutes. 
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Figure 165. Comparison of dimensionless burning rate using derived properties for 
species and orientation (1DFL8). 
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Figure 166. Comparison of burning rate using derived properties for species and 
orientation (1DFL8). 
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Test (1DFL9), Douglas Fir, along grain at 35 kW/m” for 25 minutes. 
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Figure 167. Rate of heat release and heat of combustion (1DFL9). 
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Figure 168. Comparison of dimensionless burning rate using derived properties for 
species and orientation (1DFL9). 
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Figure 169. Comparison of burning rate using derived properties for species and 
orientation (1DFL9). 
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Test (1DFX1), Douglas Fir, across grain at 75 kW/m” for 25 minutes. 
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Figure 170. Rate of heat release and heat of combustion (1DFX1). 


900 
800 4mm 
700 
600 12 mm 
ey 
= 
e 500 
he 
s 
= 
2 400 
£ 
.-*) 
f= 
300 | 
24 mm 
200% 50 mm 
ne 36mm 
100 


age grant Ambient 


0 120 240 360 480 600 720 840 960 1080 1200 1320 1440 81560 


Figure 171. Temperatures measured in sample (1DFX1). 
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Figure 172. Comparison of dimensionless burning rate using derived properties for 
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Figure 173. Comparison of burning rate using derived properties for species and 


orientation (1DFX1). 
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Figure 174. Comparison of measured and calculated thermal penetration depth (1DFX1). 


Depth [m] 


0.030 
| © Thermocouple data at Tig 
| x Thermocouple data at Tchar 

0.025 — —— Experimental data (1) 
— — — Short-time solution (2) 
— —Long-time solution (3) 

0.020 

0.015 

0.010 4 

0.005 4 

0.000 


0 120 240 360 480 600 720 840 960 1,080 1,200 1,320 1,440 1,560 - 


Figure 175. Comparison of estimated char depth from experimental data and the 


calculated char depth (1DFX1). 
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Test (1DFX2), Douglas Fir, across grain at 75 kW/m’ for 25 minutes. 
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Figure 176. Rate of heat release and heat of combustion (1DFX2). 
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Figure 177. Temperatures measured in sample (1 DFX2). 
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Figure 178. Comparison of dimensionless burning rate using derived properties for 
species and orientation (1DFX2). 
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Figure 179. Comparison of burning rate using derived properties for species and 
orientation (1DFX2). 
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Figure 180. Comparison of measured and calculated thermal penetration depth (1DFX2). 
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Figure 181. Comparison of estimated char depth from experimental data and the 
calculated char depth (1 DFX2). 
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Test (1DFX3), Douglas Fir, across grain at 75 kW/m’ for 75 minutes. 
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Figure 182. Rate of heat release and heat of combustion (1DFX3). 
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Figure 183. Temperatures measured in sample (1DFX3). 
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Figure 184. Comparison of dimensionless burning rate using derived properties for 
species and orientation (1DFX3). 
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Figure 185. Comparison of burning rate using derived properties for species and 
orientation (1 DFX3). | 
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Figure 186. Comparison of measured and calculated thermal penetration depth (1 DFX3). 
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Figure 187. Comparison of estimated char depth from experimental data and the 


calculated char depth (1DFX3). 
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Test (1DFX4), Douglas Fir, across grain at 25 kW/m? for 25 minutes. 
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Figure 188. Rate of heat release and heat of combustion (1DFX4). 
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Figure 189. Temperatures measured in sample (1DFX4). 
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Figure 190. Comparison of dimensionless burning rate using derived properties for 
species and orientation (1DFX4). 
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Figure 191. Comparison of burning rate using derived properties for species and 
orientation (1DFX4). 
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Figure 192. Comparison of measured and calculated thermal penetration depth (1DFX4). 
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Figure 193. Comparison of estimated char depth from experimental data and the 
calculated char depth (1DFX4). | 
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Test (IDFX5), Douglas Fir, across grain at 75 kW/m’ for 25 minutes. 
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Figure 194. Rate of heat release and heat of combustion (1DFX5). 
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Figure 195. Temperatures measured in sample (1DFX5). 
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Figure 196. Comparison of dimensionless burning rate using derived properties for 
species and orientation (1DFX5). 
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Figure 197. Comparison of burning rate using derived properties for species and 
orientation (1 DFX5). 
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Figure 199. Comparison of estimated char depth from experimental data and the 
calculated char depth (1DFX5). 
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Test (1DFX6), Douglas Fir, across grain at 50 kW/m? for 25 minutes. 
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Figure 200. Rate of heat release (1DFX6). 
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Figure 201. Temperatures measured in sample (1DFX6). 
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Figure 202. Comparison of dimensionless burning rate using derived properties for 
species and orientation (1DFX6). 
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Figure 203. Comparison of burning rate using derived properties for species and 
orientation (1DFX6). 
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Figure 204. Comparison of measured and calculated thermal penetration depth (1DFX6). 
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Figure 205. Comparison of estimated char depth from experimental data and the 


calculated char depth (1DFX6). 
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Test (1DFX7), Douglas Fir, across grain at 50 kW/m’ for 25 minutes. 
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Figure 206. Rate of heat release (1 DFX7). 
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Figure 207. Temperatures measured in sample (1DFX7). 
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Figure 208. Comparison of dimensionless burning rate using derived properties for 
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species and orientation (1DFX7). 
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Figure 209. Comparison of burning rate using derived properties for species and 


orientation (1 DFX7). 
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Figure 210. Comparison of measured and calculated thermal penetration depth (1 DFX7). 
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Figure 211. Comparison of estimated char depth from experimental data and the 


‘calculated char depth (1DFX7). 
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Test (1RL1), Redwood, along grain at 75 kW/m’ for 25 minutes. 
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Figure 212. Rate of heat release and heat of combustion (1RL1). 
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Figure 213. Temperatures measured in sample (1RL1). 
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Figure 214. Comparison of dimensionless burning rate using derived properties for 
species and orientation (1RL1). 


0.025 , , (2) 
Experimenta! data (1) 
- — — Small-time solution (2) 
| we 2 
0.020 J »! —  —Long-time solution (3) 
i 
i 
1 
ri ot 
é 0.015 \\ 
2 \ } (1) 
e 
re \ 
3 \ 
= 0.010 
5 \ 
a 
N 
\ 
0.005 + ict 
0S ns ee Oe a @) 
0.000 ! 
0 60 120 180 240 300 360 420 480 540 600 


Figure 215. Comparison of burning rate using derived properties for species and 
orientation (1RL1). 
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Figure 216. Comparison of measured and calculated thermal penetration depth (1RL1). 


0.030 fs 


Q) 


0.025 


0.020 


0.015 


Depth [m] 


0.010 


—— Experimental data (1) 
— — — Small-time solution (2) 
Long-time solution (3) 

© Thermocouple data at Tig 
x Thermocouple data at Tchar 


0 120 240 360 480 600 720 840 960 1,080 1,200 1,320 1,440 1,560 


Figure 217. Comparison of estimated char depth from experimental data and the 
calculated char depth (1RL1). 
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Test (1RL2), Redwood, along grain at 75 kW/m” for 25 minutes. 
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Figure 218. Rate of heat release and heat of combustion (1RL2). 
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Figure 219. Temperatures measured in sample (1RL2). 
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Figure 220. Comparison of dimensionless burning rate using derived properties for 
species and orientation (1RL2). 
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Figure 221. Comparison of burning rate using derived properties for species and 
orientation (1RL2). 


276 


1.0 


Depth [m] 


| 
0.045 
0.040 
0.035 Oo 
0.030 | D eatin 
0.025 | os Eee 
¢ ‘Spies os 


Pre-ignition solution (1) 
— — — Small-time solution (2) 
— —Long-time solution (3) 
tO Thermocouple data 


0 60 120 180 240 300 360 ; 420 


Time |s] 


Figure 222. Comparison of measured and calculated thermal penetration depth (1RL2). 
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Figure 223. Comparison of estimated char depth from experimental data and the 


calculated char depth (1RL2). 
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Test (1RL3), Redwood, along grain at 25 kW/m” for 25 minutes. 
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Figure 224. Rate of heat release and heat of combustion (1RL3). 
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Figure 225. Temperatures measured in sample (1RL3). 
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Figure 226. Comparison of dimensionless burning rate using derived properties for 
species and orientation (1RL3). 
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Figure 227. Comparison of burning rate using derived properties for species and 


orientation (1RL3). 
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Figure 228. Comparison of measured and calculated thermal penetration depth (1RL3). 
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Figure 229. Comparison of estimated char depth from experimental data and the 
calculated char depth (1RL3). 
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Test (1RLA4), Redwood, along grain at 75 kW/m” for 25 minutes. 
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Figure 230. Rate of heat release and heat of combustion (1RL4). 
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Figure 231. Temperatures measured in sample (1RL4). 
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Experimental data (1) 
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Figure 232. Comparison of dimensionless burning rate using derived properties for 
species and orientation (1RL4). 
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Figure 233. Comparison of burning rate using derived properties for species and 
orientation (1RL4). 
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Figure 234. Comparison of measured and calculated thermal penetration depth (1RL4). 
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Figure 235. Comparison of estimated char depth from experimental data and the 


calculated char depth (1RL4). 
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Test (1RL6), Redwood, along grain at 50 kW/m’ for 25 minutes. 
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Figure 236. Rate of heat release (1RL6). 
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Figure 237. Temperatures measured in sample (1RL6). 
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Figure 238. Comparison of dimensionless burning rate using derived properties for 
species and orientation (1RL6). 
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Figure 239. Comparison of burning rate using derived properties for species and 
orientation (1RL6). 
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Figure 240. Comparison of measured and calculated thermal penetration depth (1RL6). 
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Figure 241. Comparison of estimated char depth from experimental data and the 


calculated char depth (1RL6). 
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Test (1RL7), Redwood, along grain at 50 kW/m? for 25 minutes. 
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Figure 242. Rate of heat release (1RL7). 
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Figure 243. Temperatures measured in sample (1RL7). 
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Figure 244. Comparison of dimensionless burning rate using derived properties for 
species and orientation (1RL7). 
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Figure 245. Comparison of burning rate using derived properties for species and 
orientation (1RL7). 
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Figure 246. Comparison of measured and calculated thermal penetration depth (1RL7). 
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Figure 247. Comparison of estimated char depth from experimental data and the 


calculated char depth (1RL7). 
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Test (1RL8), Redwood, along grain at 50 kW/m’ for 25 minutes. 
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Figure 248. Rate of heat release and heat of combustion (1RL8). 
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Figure 249. Comparison of dimensionless burning rate using derived properties for 
species and orientation (1RL8). 
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Figure 250. Comparison of burning rate using derived properties for species and 
orientation (1RL8). 
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Test (IRL9), Redwood, along grain at 35 kW/m” for 25 minutes. 
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Figure 251. Rate of heat release and heat of combustion (1RL9). 
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Figure 252. Comparison of dimensionless burning rate using derived properties for 
species and orientation (1RL9). 
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Figure 253. Comparison of burning rate using derived properties for species and 
orientation (1RL9). 
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Test (1RX1), Redwood, across grain at 75 kW/m’? for 25 minutes. 
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Figure 254. Rate of heat release and heat of combustion (1RX1). 
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Figure 255. Temperatures measured in sample (1RX1). 
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Figure 256. Comparison of dimensionless burning rate using derived properties for 
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Figure 257. Comparison of burning rate using derived properties for species and 


orientation (1RX1). 


295 


Depth [m] 


Pre-ignition solution (1) 
— — - Short-time solution (2) 
— —Long-time solution (3) 
0). «Thermocouple data 


0 30 60 90 120 150 180 210 240 270 300 330 
Time [s] 


Figure 258. Comparison of measured and calculated thermal penetration depth (1RX1). 
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Figure 259. Comparison of estimated char depth from experimental data and the 


calculated char depth (1RX1). 
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Test (1RX2), Redwood, across grain at 75 kW/m’ for 25 minutes. 
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Figure 260. Rate of heat release and heat of combustion (1RX2). 
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Figure 261. Temperatures measured in sample (1RX2). 
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Experimental data (1) 
— — — Small-time solution (2) 
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Figure 262. Comparison of dimensionless burning rate using derived properties for 
species and orientation (1RX2). 
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Figure 263. Comparison of burning rate using derived properties for species and 
orientation (1RX2). 
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Figure 264. Comparison of measured and calculated thermal penetration depth (1RX2). 
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Figure 265. Comparison of estimated char depth from experimental data and the 


calculated char depth (1RX2). 
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Test (IRX3), Redwood, across grain at 25 kW/m” for 25 minutes. 
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Figure 266. Rate of heat release and heat of combustion (1RX3). 
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Figure 267. Temperatures measured in sample (1RX3). 
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Figure 268. Comparison of dimensionless burning rate using derived properties for 


species and orientation (1RX3). 
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Figure 269. Comparison of burning rate using derived properties for species and 
orientation (1RX3). 
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Figure 270. Comparison of measured and calculated thermal penetration depth (1RX3). 
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Figure 271. Comparison of estimated char depth from experimental data and the 
calculated char depth (1RX3). 
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Test (1RX4), Redwood, across grain at 25 kW/m? for 25 minutes. 
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Figure 272. Rate of heat release and heat of combustion (1RX4). 
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Figure 273. Temperatures measured in sample (1RX4). 
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Figure 274. Comparison of dimensionless burning rate using derived properties for 
species and orientation (1RX4). 
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Figure 275. Comparison of burning rate using derived properties for species and 
orientation (1RX4). 
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Figure 276. Comparison of measured and calculated thermal penetration depth (1RX4). 
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Figure 277. Comparison of estimated char depth from experimental data and the 


calculated char depth (1RX4). 
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Test (1RX6), Redwood, across grain at 25 kW/m? for 75 minutes. 
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Figure 278. Rate of heat release (1RX6). 
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Figure 279. Temperatures measured in sample (1RX6). 
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Figure 280. Comparison of dimensionless burning rate using derived properties for 
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Figure 281. Comparison of burning rate using derived properties for species and 


orientation (1RX6). 
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Figure 282. Comparison of measured and calculated thermal penetration depth (1RX6). 
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Figure 283. Comparison of estimated char depth from experimental data and the 


calculated char depth (1RX6). 
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Test (1RX7), Redwood, across grain at 50 kW/m’ for 25 minutes. 
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Figure 284. Rate of heat release (1RX7). 
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Figure 285. Temperatures measured in sample (1RX7). 
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Figure 286. Comparison of dimensionless burning rate using derived properties for 
species and orientation (1RX7). 
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Figure 287. Comparison of burning rate using derived properties for species and 
orientation (1RX7). 
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Figure 288. Comparison of measured and calculated thermal penetration depth (1RX7). 
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Figure 289. Comparison of estimated char depth from experimental data and the 


calculated char depth (1RX7). 


311 


Test (IRX8), Redwood, across grain at 50 kW/m? for 25 minutes. 
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Figure 290. Rate of heat release (1RX8). 
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Figure 291. Temperatures measured in sample (1RX8). 
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Figure 292. Comparison of dimensionless burning rate using derived properties for 


species and orientation (1RX8). 
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Figure 293. Comparison of burning rate using derived properties for species and 
orientation (1RX8). 
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Figure 294. Comparison of measured and calculated thermal penetration depth (1RX8). 
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Figure 295. Comparison of estimated char depth from experimental data and the 


calculated char depth (1RX8). 
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Test (1RX9), Redwood, across grain at 50 kW/m’ for 25 minutes. 
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Figure 296. Rate of heat release and heat of combustion (1RX9). 
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Figure 297. Comparison of dimensionless burning rate using derived properties for 
species and orientation (1RX9). 
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Figure 298. Comparison of burning rate using derived properties for species and 
orientation (1RX9). 
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Test (1RX10), Redwood, across grain at 35 kW/m’ for 25 minutes. 
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Figure 300. Comparison of dimensionless burning rate using derived properties for 
species and orientation (1RX10). 
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Figure 301. Comparison of burning rate using derived properties for species and 
orientation (1RX10). 
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Test (10L1), Red oak, along grain at 25 kW/m? for 25 minutes. 
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Figure 302. Rate of heat release and heat of combustion (10L1). 
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Figure 303. Temperatures measured in sample (10L1). 
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Figure 304. Comparison of dimensionless burning rate using derived properties for 
species and orientation (1OL1). 
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Figure 305. Comparison of burning rate using derived properties for species and 
orientation (1OL1). 
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Figure 306. Comparison of measured and calculated thermal penetration depth (10L1). 
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Figure 307. Comparison of estimated char depth from experimental data and the 


calculated char depth (10OL1). 
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Test (10L2), Red oak, along grain at 25 kW/m? for 75 minutes. 
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Figure 308. Rate of heat release (10L2). 
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Figure 309. Temperatures measured in sample (10L2). 
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Figure 310. Comparison of dimensionless burning rate using derived properties for 


species and orientation (10L2). 
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Figure 311. Comparison of burning rate using derived properties for species and 
orientation (1OL2). 
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Figure 313. Comparison of estimated char depth from experimental data and the 
calculated char depth (10L2). 
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Test (10L3), Red oak, along grain at 75 kW/m’ for 25 minutes. 
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Figure 315. Temperatures measured in sample (10L3). 
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Figure 316. Comparison of dimensionless burning rate using derived properties for 
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species and orientation (10L3). 
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Figure 317. Comparison of burning rate using derived properties for species and 


orientation (10L3). 
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Figure 318. Comparison of measured and calculated thermal penetration depth (10L3). 
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Figure 319. Comparison of estimated char depth from experimental data and the 


calculated char depth (10L3). 
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Figure 320. Rate of heat release and heat of combustion (10L4). 
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Figure 321. Temperatures measured in sample (10L4). 
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Figure 322. Comparison of dimensionless burning rate using derived properties for 
species and orientation (10L4). 
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Figure 323. Comparison of burning rate using derived properties for species and 
orientation (10L4). 
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Figure 324. Comparison of measured and calculated thermal penetration depth (10L4). 
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Figure 325. Comparison of estimated char depth from experimental data and the 


calculated char depth (10L4). 
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Test (10L5), Red oak, along grain at 75 kW/m’ for 25 minutes. 
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Figure 326. Rate of heat release and heat of combustion (10L5). 
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Figure 327. Temperatures measured in sample (10L5). 
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Figure 328. Comparison of dimensionless burning rate using derived properties for 


0.020 


0.015 


Burning rate [kg/s.m’] 


0.010 


0.005 


species and orientation (10L5). 


Expenmental data (1) 
— — — Small-time solution (2) 
— —Long-time solution (3) 


60 120 180 240 300 360 420 480 540 600 
Time [s] 


Figure 329. Comparison of burning rate using derived properties for species and 


orientation (1OL5). 
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Figure 330. Comparison of measured and calculated thermal penetration depth (10L5). 
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Figure 331. Comparison of estimated char depth from experimental data and the 


calculated char depth (1OL5S). 
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Figure 332. Rate of heat release (10X1). 
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Figure 333. Temperatures measured in sample (10X1). 
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Figure 334. Comparison of dimensionless burning rate using derived properties for 
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Figure 335. Comparison of burning rate using derived properties for species and 


orientation (10X1). 


335 


Pre-ignition solution (1) 
~ — — Short-time solution (2) 
— —Long-time solution (3) 


Depth (m] 


0 30 60 90 120 150 180 210 240 270 300 330 360 390 420 


Figure 336. Comparison of measured and calculated thermal penetration depth (10X1). 
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Figure 337. Comparison of estimated char depth from experimental data and the 
‘calculated char depth (10X1). 
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Test (10X2), Red oak, along grain at 25 kW/m” for 75 minutes. 
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Figure 338. Rate of heat release (10X2). 
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Figure 339. Temperatures measured in sample (10X2). 
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Figure 340. Comparison of dimensionless burning rate using derived properties for 
species and orientation (10X2). 
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Figure 341. Comparison of burning rate using derived properties for species and 
orientation (1OX2). 
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Figure 342. Comparison of measured and calculated thermal penetration depth (10X2). 
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Figure 343. Comparison of estimated char depth from experimental data and the 


calculated char depth (10X2). 
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Test (10X3), Red oak, along grain at 75 kW/m’ for 25 minutes. 
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Figure 344. Rate of heat release and heat of combustion (10X3). 
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Figure 345. Temperatures measured in sample (10X3). 
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Figure 346. Comparison of dimensionless burning rate using derived properties for 
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Figure 347. Comparison of burning rate using derived properties for species and 


orientation (10X3). 
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Figure 348. Comparison of measured and calculated thermal penetration depth (10X3). 
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Figure 349. Comparison of estimated char depth from experimental data and the 


calculated char depth (10X3). 
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Test (10X4) : Red oak, along grain at 75 kW/m’ for 25 minutes. 
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Figure 350. Rate of heat release and heat of combustion (10X4). 
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Figure 351. Temperatures measured in sample (10X4). 
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Figure 352. Comparison of dimensionless burning rate using derived properties for 
species and orientation (10OX4). 
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Figure 353. Comparison of burning rate using derived properties for species and 
orientation (10X4). 
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Figure 354. Comparison of measured and calculated thermal penetration depth (10X4). 
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Figure 355. Comparison of estimated char depth from experimental data and the 
calculated char depth (10X4). 
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Test 10X5 : Red oak, along grain at 75 kW/m? for 25 minutes. 


' 30 


asa Rate of heat release 
—— Heat of combustion 25 


Rate of heat release [kW/m 7] 
Heat of combustion [MJ/kg] 


360 480 600 720 840 960 1,080 1,200 1,320 1,440 
Time [s] 


Figure 356. Rate of heat release and heat of combustion (10X5). 


800 


4mm 


12 mm 


Temperature [°C] 


36 mm 
50 mm 


Ambient 
0 ee 
0 120 240 360 480 600 720 840 960 1080 1200 1320 1440 #1560 


Figure 357. Temperatures measured in sample (10X5). 
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Figure 358. Comparison of dimensionless burning rate using derived properties for 
species and orientation (10X5). 
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Figure 359. Comparison of burning rate using derived properties for species and 
orientation (10OX5). 
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Figure 360. Comparison of measured and calculated thermal penetration depth (10X5). 
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Figure 361. Comparison of estimated char depth from experimental data and the 
calculated char depth (1OXS5). 
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Test (1ML1), Maple, along grain at 75 kW/m? for 25 minutes. 
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Figure 362. Rate of heat release and heat of combustion (1ML1). 
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Figure 363. Temperatures measured in sample (1ML1). 
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Figure 364. Comparison of dimensionless burning rate using derived properties for 
species and orientation (1ML1). 
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Figure 365. Comparison of burning rate using derived properties for species and 
orientation (1ML1). 
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Figure 366. Comparison of measured and calculated thermal penetration depth (1ML1). 
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Figure 367. Comparison of estimated char depth from experimental data and the 


calculated char depth (1ML1). 
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Test (1ML2), Maple, along grain at 75 kW/m? for 25 minutes. 
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Figure 368. Rate of heat release and heat of combustion (1ML2). 
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Figure 369. Temperatures measured in sample (1ML2). 
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Figure 370. Comparison of dimensionless burning rate using derived properties for 
species and orientation (1ML2). 
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Figure 371. Comparison of burning rate using derived properties for species and 
orientation (1ML2). 
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Figure 372. Comparison of measured and calculated thermal penetration depth (1ML2). 
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Figure 373. Comparison of estimated char depth from experimental data and the 
calculated char depth (1ML2). 
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Test (1ML3), Maple, along grain at 75 kW/m? for 25 minutes. 
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Figure 374. Rate of heat release and heat of combustion (1ML3). 
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Figure 375. Temperatures measured in sample (1ML3). 
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Figure 376. Comparison of dimensionless burning rate using derived properties for 
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Figure 377. Comparison of burning rate using derived properties for species and 


orientation (1ML3). 
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Figure 378. Comparison of measured and calculated thermal penetration depth (1ML3). 
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Figure 379. Comparison of estimated char depth from experimental data and the 


calculated char depth (1ML3). 
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Test (1ML4), Maple, along grain at 25 kW/m’ for 25 minutes. 
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Figure 380. Rate of heat release and heat of combustion (1ML4). 
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Figure 381. Temperatures measured in sample (1ML4). 
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Figure 382. Comparison of dimensionless burning rate using derived properties for 
species and orientation (1ML4). 
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Figure 383. Comparison of burning rate using derived properties for species and 
orientation (1ML4). 
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Figure 384. Comparison of measured and calculated thermal penetration depth (1ML4). 
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Figure 385. Comparison of estimated char depth from experimental data and the 


calculated char depth (1ML4). 
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Test (1ML5), Maple, along grain at 25 kW/m? for 75 minutes. 
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Figure 386. Rate of heat release (1ML5S). 
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Figure 387. Temperatures measured in sample (1ML5). 


361 


Test (1MX1), Maple, across grain at 25 kW/m? for 25 minutes. 
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Figure 388. Rate of heat release (1MX1). 
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Figure 389. Temperatures measured in sample (1MX1). 
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Figure 390. Comparison of dimensionless burning rate using derived properties for 
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Figure 391. Comparison of burning rate using derived properties for species and 


orientation (1MX1). 
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Figure 392. Comparison of measured and calculated thermal penetration depth (1MX1). 
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Figure 393. Comparison of estimated char depth from experimental data and the 


calculated char depth (1MX1). 
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Test (1MX2), Maple, across grain at 25 kW/m” for 75 minutes. 
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Figure 394. Rate of heat release (1MX2). 
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Figure 395. Temperatures measured in sample (1MX2). 
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Figure 396. Comparison of dimensionless burning rate using derived properties for 
species and orientation (1MX2). 
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Figure 397. Comparison of burning rate using derived properties for species and 
orientation (1MX2). 
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Figure 398. Comparison of measured and calculated thermal penetration depth (1MX2). 
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Figure 399. Comparison of estimated char depth from experimental data and the 


calculated char depth (1MX2). 
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Test (1MX3), Maple, across grain at 75 kW/m” for 25 minutes. 
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Figure 400. Rate of heat release and heat of combustion (1MX3). 
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Figure 401. Temperatures measured in sample (1MX3). 
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Figure 402. Comparison of dimensionless burning rate using derived properties for 
species and orientation (1 MX3). 
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Figure 403. Comparison of burning rate using derived properties for species and 
orientation (1MX3). 
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Figure 404. Comparison of measured and calculated thermal penetration depth (1 MX3). 
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Figure 405. Comparison of estimated char depth from experimental data and the 


calculated char depth (1MX3). 
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Test (IMX4), Maple, across grain at 75 kW/m’ for 25 minutes. 
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Figure 406. Rate of heat release and heat of combustion (1MX4). 
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Figure 407. Temperatures measured in sample (1MX4). 
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Figure 408. Comparison of dimensionless burning rate using derived properties for 
species and orientation (1MX4). 
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Figure 409. Comparison of burning rate using derived properties for species and 
orientation (1MX4). 
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Figure 410. Comparison of measured and calculated thermal penetration depth (1MX4). 
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Figure 411. Comparison of estimated char depth from experimental data and the 
calculated char depth (1MX4). 
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Test (1MXS), Maple, across grain at 25 kW/m’? for 75 minutes. 
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Figure 412. Rate of heat release and heat of combustion (1MXS). 
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Figure 413. Temperatures measured in sample (1MX5). 
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Figure 414. Comparison of dimensionless burning rate using derived properties for 
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Figure 415. Comparison of burning rate using derived properties for species and 


orientation (1MXS5). 
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Figure 416. Comparison of measured and calculated thermal penetration depth (1MX5). 
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Figure 417. Comparison of estimated char depth from Deters! data and the 
calculated char depth (1MXS5). 
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